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ABSTRACT 

This  thesis  devoted  to  studying  the  properties  of  a  digital  M-ary  signaling  scheme  called 
MFSK'QPSK  which  combines  the  features  of  M-ary  FSK  and  QPSK  modulation.  This 
scheme  is  shown  to  be  a  form  of  biorthogonal  signaling,  and  a  receiver  that  is  optimum 
in  error  rate  sense  is  introduced  and  analyzed  in  order  to  obtain  its  bit  error  rate  (BER) 
performance  when  signal  interference  can  be  modeled  as  additive  white  Gaussian  noise. 
The  bandwidth  of  the  M  FSK;  QPSK  signaling  scheme  is  derived  and  compared  with  that 
of  M-ar>'  PSK  and  M-ar\'  FSK  signaling.  MFSK/QPSK  signaling  results  in  receivers 
that  provide  error  rates  much  lower  than  corresponding  receivers  for  M-ary  PSK  and 
slightly  lower  than  those  for  M-ar>'  FSK,  while  maintaining  a  bandwidth  efficiency  at 
least  three  times  greater  than  that  of  M-ar\-  FSK  though  not  as  great  as  that  of  M-ary 
PSK  signaling.  A  significant  effort  has  been  devoted  to  examining  a  detection  method 
known  as  direct  bit  detection  (DBD)  which  has  been  applied  to  M-arv'  PSK  and 
M FSK  QPSK  signaling.  DBD  receivers  are  introduced,  and  the  BER  of  the 
M  FSK  QPSK  DBD  receiver  is  derived  and  compared  with  that  of  the  optimum  (symbol 
discrimination)  receiver.  Lastly,  the  efiects  on  performance  due  to  Rayleigh  fading  are 
analyzed  and  evaluated  for  M-ar>"  PSK  DBD  receivers  as  well  as  for  the  MFSK  QPSK 
optimum  (symbol  discrimination)  and  DBD  receiver. 
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I.     INTRODUCTION 

A.  BACKGROUND 

There  exists  an  increasing  need  for  high-data-rate,  reliable  digital  communication 
systems.  The  reasons  for  this  need  stem  from  the  desire  by  an  increased  number  of  users 
to  transmit  larger  amounts  of  data  via  digital  communication  schemes  which  provide 
more  reliability,  flexibility,  and  lower  cost  when  compared  with  analog  communication 
systems.  These  increased  demands  come  at  a  time  when  the  available  bandwidth  in  the 
free  space  channel  is  becoming  a  limited  resource.  These  data  transmission  requirements 
coupled  with  limited  channel  resources  are  the  driving  force  behind  the  continual  search 
for  digital  modulation  techniques  which  operate  at  low  error  rates  and  are  bandwidth 
efficient.  Present  signaling  schemes  in  use  are  the  well  known  M-ar\'  Phase  Shift  Keying 
(M-PSK)  and  M-ar>'  Frequency  Shift  Keying  (M-FSK).  These  signaling  schemes  are 
briefly  described  in  Sections  C  and  D  of  this  chapter,  and  the  corresponding  receiver 
structures  for  recovering  the  digital  information  are  described  in  Section  E,  while  the 
performance  of  these  receivers,  specified  in  terms  of  the  symbol  error  rates,  for  both 
M-ar\'  PSK  and  M-an.'  FSK  modulation  are  given  in  Section  F.  The  bandwidth  re- 
quirements associated  with  these  two  modulation  schemes  are  given  in  Section  G. 

B.  CONTENTS  AND  SUMMARY  OF  RESULTS 

This  thesis  explores  a  hybrid  form  of  digital  modulation,  namely  MFSK'QPSK, 
which  embodies  features  of  both  M-ar\-  FSK  and  Quaternan,'  (or  four  phase)  PSK.  The 
analysis  compares  and  contrasts  the  properties  of  this  modulation  scheme  to  the  widely 
known  M-ar\'  PSK  and  M-ar\-  FSK  modulation.  In  Chapter  3,  MFSK,QPSK  is  shown 
to  be  a  form  of  biorthogonal  signaling,  and  therefore  its  corresponding  receiver  per- 
formance is  easily  revealed  given  what  is  already  known  for  this  class  of  signals.  The 
optimum  MFSK'QPSK  receiver  that  recovers  the  digital  information  with  minimum 
probability  of  error  under  specific  noise  interference  assumptions  is  analyzed  from  a 
performance  standpoint,  which  traditionally  has  been  given  in  terms  of  the  Symbol  Error 
Rate  (SER).  Because  a  more  meaningful  performance  parameter  is  the  receiver's  Bit 
Error  Rate  (BER),  Chapter  3  has  been  devoted  to  deriving  the  receiver's  BER  as  a 
function  of  signal-to-noise  ratio,  along  with  the  bandwidth  requirements  and  the  band- 
width efficiency  of  the  MFSK/QPSK  scheme.  Comparisons  are  made  between  the 
MFSK/QPSK  scheme  and  conventional  phase  and  frequency  modulation  methods  with 


respect  to  SER  performance  and  bandwidth  efficiency  in  Chapter  5.  The  results  show 
that  the  MFSK/QPSK  scheme  when  compared  to  M-ary-  PSK,  does  not  attain  the  high 
bandwidth  efficiencies  achieved  by  M-ar\'  PSK,  yet  outperforms  M-ary  PSK  substan- 
tially in  terms  of  SER  performance,  while  performing  slightly  better  than  M-ary  FSK  in 
terms  of  SER  performance  and  only  requiring  one  third  to  one  fourth  the  bandwidth 
required  by  an  M-ar>'  FSK  system. 

Additionally,  this  thesis  addresses  the  properties  of  a  class  of  digital  communication 
receivers  which  decode  the  individual  bits  that  make  up  a  symbol  directly,  without  need 
for  symbol-to-bit  mapping  post-detection  circuitr\'.  These  are  referred  to  as  direct  bit 
detection  (DBD)  receivers,  because  each  of  the  bits  that  make  up  a  symbol  are  deter- 
mined independently.  These  receivers  have  been  proposed  and  analyzed  for  M-ary  PSK 
Ref  1  [pp.  461-463]  as  well  as  for  the  MFSK/QPSK  signaling  scheme  Ref  2  [pp.  16-101] 
under  discussion  here.  DBD  receiver  structures  for  M-ary  PSK  are  introduced  in 
Chapter  2,  and  previously  derived  expressions  for  their  BER  performance  are  presented. 
In  Chapter  3,  a  DBD  receiver  structure  for  MFSK/QPSK  signaling  is  described,  and 
analysis  is  carried  out  in  order  to  derive  its  BER  performance,  with  corresponding  results 
presented  in  graphical  form  in  Chapter  5.  Also  in  Chapter  5  the  BER  performance  of 
the  MFSK/QPSK  DBD  receiver  is  compared  to  that  of  the  optimum  receiver  that  per- 
forms symbol  detection  followed  by  symbol-to-bit  mappings  prior  to  delivering  binary 
data  to  the  intended  user.  Furthermore,  the  results  demonstrate  that  while  the  DBD 
receiver  does  not  perform  as  well  as  the  optimum  receiver  in  BER  sense,  its  simplified 
receiver  structure  and  direct  bit  detection  properties  offer  other  potential  advantages. 

In  choosing  a  modulation  method,  it  is  important  that  the  effect  of  system  disturb- 
ances such  as  Rayleigh  fading,  which  is  common  in  digital  microwave  radio  transmission 
applications,  be  analyzed  in  order  to  fully  understand  all  aspects  of  expected  system 
performance  for  the  different  signaling  schemes  under  consideration.  The  effect  of 
Rayleigh  fading  on  the  BER  performance  of  the  MFSK/QPSK  receivers  as  well  as  of  the 
M-ar\^  PSK  DBD  receivers  is  analyzed  and  evaluated  in  Chapter  4.  The  results  demon- 
strate that  in  all  cases,  Rayleigh  fading  can  severely  degrade  receiver  performance. 

C.     MARY  PHASE-SHIFT  KEYING 

M-ary  signaling  techniques  involve  forming  groups  of  k  bits  to  generate  iV/=2* 
symbols  which  are  assigned  to  M  different  waveforms  suitable  for  transmission  over  a 
channel.  One  form  of  M-ar\-  signaling,  namely  M-ary  Phase-Shift  Keying  (M-PSK),  is 
a   bandwidth   efficient    digital   modulation   scheme.      The   phase   of  the   transmitted 


sinusoidal  carrier  waveform  can  take  on  one  of  M  different  phases  to  represent  M  dif- 
ferent svmbols.   The  M-arv'  PSK  waveform  can  be  described  mathematically  as 


s,{i)=J^  cos[2;:/o;  +  ^]  (1.1) 

for  i  =  0,  I,.",  M-1  and  0  <  t  <T„  where  generally  (as  is  the  case  here)  the  phase  dif- 
ference between  waveform  pairs  is  taken  to  be  constant.  Each  symbol  or  waveform 
Si{t)  is  transmitted  at  a  rate  of  R,  =  -jr  symbols/sec,  where  T,  is  the  time  duration  of  each 
waveform  sit).  For  k  =  2  and  A/=4,  the  resulting  modulation  form  is  known  as 
Quaternar\'  Phase-Shift  Keying  (QPSK),  where  the  phase  difference  between  waveform 
pairs  is  —  radians.  For  8-PSK,  in  which  k  =  3  and  M  =  8,  waveform  pairs  are  separated 
by  a  phase  of—  radians.  Signal  constellations  for  QPSK  and  8-PSK  as  well  as  symbol 
assignments  for  each  signal  (or  phasor  representation)  are  presented  in  Figure  1  on  page 
4. 

D.    M-ARY  FREQUENCY-SHIFT  KEYING 

M-ar\'  Frequency-Shift  Keying  (M-FSK)  is  another  form  of  M-ar\'  signaling,  in 
which  the  waveforms  are  mathematically  expressed  as 


5,(0=^^  coslnf^t  (1.2) 

where  the  frequencies/,  typically  are  given  by 

y;=/o  +  /A/  (1.3) 

for  i  =  0,  1,...,  M-1  and  0  <  t  <  T,.  The  signal  duration  is  normally  chosen  so  that/r, 
is  an  integer  or  at  least /r,  >  >  1  holds.  The  frequency  separation  A/ is  often  set  equal 
to  R,  providing  minimum  frequency  spacing  between  signals  s,{t)  ,  while  forming  a  set 
of  orthogonal  signals.  For  /:  =  2  and  M  =  A,  the  resulting  modulation  is  Quadrature 
Frequency- Shift  Keying  (QFSK).  Magnitude  spectra  for  QFSK  and  M-ar\'  FSK  (  for 
arbitrary  M  )  are  shown  in  Figure  2  on  page  5  and  in  Figure  3  on  page  6,  respectively. 

E.     RECEIVERS  OF  M-ARY  PSK  AND  M-ARY  FSK 

Signals  transmitted  through  a  physical  channel  always  encounter  interference.  The 
SER  results  given  later  in  this  chapter  and  the  BER  analysis  presented  in  Chapter  3  are 
obtained  under  the  assumption  that  this  interference  can  be  modeled  as  a  sample  func- 
tion of  an  additive  white  Gaussian  noise  (AWGN)  process  having  double  sided  power 
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Figure  1.       Signal  Constellations  for  a)  QPSK  and  b)  8-PSK 
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spectral  density  (PSD)  of  — r-.  Lndcr  this  assumption,  receivers  can  be  derived  that  are 
optimum  in  the  sense  that  the  probabihty  of  a  symbol  error  is  minimized.  One  structural 
form  of  an  optimum  receiver  for  M-ary  PSK  is  shown  in  Figure  4  on  page  8.  This  re- 
ceiver as  well  as  all  receivers  considered  herein  is  assumed  to  be  coherent,  meaning  that 
the  phase  of  the  sinusoidal  carrier  is  known  at  the  receiver.  The  concept  of  coherence 
applies  to  this  receiver  in  that  the  cosine  and  sine  functions  which  are  correlated  with  the 
received  signal  are  in  phase  with  the  unmodulated  carrier  sinusoid.  Because  the  receiver 
correlates  the  received  sisnal  with  a  set  of  sinusoidal  functions,  which  in  this  case  con- 
sists  of  each  of  the  possible  transmitted  signals,  it  is  referred  to  as  a  correlation  receiver. 
An  optimum  receiver  structure  for  M-ary  FSK  is  shown  in  Figure  5  on  page  9, 
which  also  corresponds  to  a  correlation  receiver.  The  correlator  with  a  local  signal  input 
which  most  closely  matches  the  incoming  signal,  will  produce  in  the  absence  of  interfer- 
ence, the  largest  output  so  that  the  receiver  can  decide  which  signal  was  sent  according 
to  which  correlator  output  is  the  largest. 

F.     SYMBOL  ERROR  RATES  FOR  MARY  PSK  AND  M-ARY  FSK 

For  coherent  demodulation  of  M-arv'  PSK  using  the  optimum  receiver  shown  in 
Figure  4,  the  SER  is  given  by  Ref  3  [pp.  204-207] 


Ps{e]=\--TT  +  Q\J¥-  sin   "" 


FT  (1-4) 


--^[^  ''"  exp(-v')[l  -2Q(>v2  cot-^)]^y 


where 


poo 

QiO^l    giOdC  (1.5) 

The  integral  on  the  right  hand  side  of  Eq.  1.4  has  no  known  closed  form,  and  therefore 

£ 

must  be  evaluated  numerically  if  plotting  oi'  P,{e]  as  a  function  of-rr-  is  desired.    For 

f  "  ^  ■''0 

high  values  of-rr-,  the  SER  can  be  tightly  upperbounded  by 


''''^''=^V"f ''":^  <'-'^ 
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Figure  4.       Receiver  Structure  for  M-ary  PSK 
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Figure  5.      Receiver  Structure  for  M-ar>'  FSK 


For  M-arv'  FSK  modulation  and  an  AWGN  interference  model,  the  optimum  co- 
herent receiver  of  Figure  5  on  page  9  achieves  a  probability  of  correctly  decoding  a 
symbol  given  by  Ref  3  [pp.  200-203] 


PsiC]  = 


1 


^wf 


g{x)dx 


M-] 


gWy 


(1.7) 


A' 


where  a  =  — r-,  and  ^(;)  is  the  zero  mean  unit  variance  Gaussian  distribution,  or 


^(C)  =  -T=-  expf  -  -J 


V 


(1.8) 


This  expression  can  be  rewritten  as 


/»oo_ 


PsiQ 


l-0{v  +  VJ- 


M-\ 


g{y)dy  =  1  -  P,{e] 


(1.9) 


Also,  the  SER  can  be  closely  upperbounded  for  high  values  of-rr-  by 

An 


PAe)<iM-mUi-^ 


(1.10) 


G.     BANDWIDTH  REQUIREMENTS 

For  M-ar>'  PSK,  the  carrier  is  switched  at  the  rate  of  R,  times  per  second.    Since 
there  are  k  bits  associated  with  each  symbol,  the  bit  rate  /?^  can  be  expressed  as 


Rb  =  kRs 


(1.11) 


The  null-to-null  bandwidth  IV  for  M-ary  PSK  is  given  by 


2R, 


W=2R,=     , 
k 


(1.12) 


For  M-ary  FSK  the  carrier  is  also  switched  at  the  symbol  rate  R,  .  Looking  at  Fig- 
ure 3  it  can  be  seen  that  the  total  null-to-null  bandwidth  left-most  main  lobe  to  right- 
most main  lobe  for  M-ar\-  FSK  can  be  written  as 
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W={M+\)R,  =  - ^--^  (1.13) 

From  Eq.  1.13  it  is  apparent  that  the  null-to-null  bandwidth  of  M-ary  FSK  increases 
with  increasing  values  of  k,  however  Eq.  1.12  shows  that  the  null-to-null  bandwidth  of 
M-ar>'  PSK  decreases  with  increasing  k.  This  property  makes  M-ar\'  PSK  a  desirable 
modulation  method  when  bandwidth  must  be  minimized.  In  the  next  chapter  M-ar>' 
PSK  receivers  of  a  different  structure  than  the  one  shown  in  Figure  4  are  introduced, 
and  their  bit  error  rates  presented. 
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II.     DIRECT  BIT  DETECTION  METHODS  FOR  M-ARY  PSK 

A.  MARY  PSK  DIRECT  BIT  DETECTION  RECEIVERS 

M-ary  PSK  modulation  as  discussed  in  Chapter  1,  Section  A,  is  a  signaling  scheme 
that  achieves  bandwidth  conservation  at  the  expense  of  signal-to-noise  ratio  required  to 
maintain  a  certain  bit  error  rate  performance  level.  The  two  forms  of  M-arv'  PSK  sig- 
naling discussed  in  this  chapter  are  8-PSK  and  16-PSK.  For  8-PSK,  k  =  3  bits/symbol 
and  A/ =8  symbols,  while  for  16-PSK,  k  =  A  bits/symbol  and  M  =  16  symbols.  The 
M-ar>'  PSK  receiver  discussed  in  Chapter  1,  Section  E  and  most  other  optimum  receiver 
structures  discussed  in  the  literature,  require  symbol-to-bit  mapping  to  be  carried  out 
after  a  decision  has  been  made  as  to  what  signal  was  transmitted.  A  receiver  structure 
which  does  not  require  this  symbol-to-bit  mapping  is  proposed  in  Ref  4  [pp.  3-5]  and  in 
Ref  1  [pp.  461-463]  for  8-PSK  and  16-PSK,  respectively.  Block,  diagrams  for  these  re- 
ceivers are  shown  in  Figure  6  on  page  13  for  8-PSK,  and  in  Figure  7  on  page  14  for 
16-PSK.  Because  each  of  the  bits  are  determined  independently  by  these  receivers  by 
appropriate  processing  of  the  correlator  outputs,  these  receivers  are  referred  to  as  direct 
bit  detection  (DBD)  receivers.  The  processes  used  by  DBD  receivers  for  recovering  bits 
make  it  possible  to  adjust  the  BER  of  the  individual  bits  within  a  symbol  by  means  of 
var\'ing  the  angles  between  the  signal  waveforms.  These  angle  adjustments  are  shown 
in  Figure  8  on  page  15  for  8-PSK,  and  in  Figure  9  on  page  16  for  16-PSK,  where  a  and 
y6  are  the  adjustable  phase  parameters. 

B.  PERFORMANCE  OF  DBD  RECEIVERS  FOR  MARY  PSK 

A  BER  analysis  for  the  8-PSK  DBD  receiver  was  carried  out  by  Ref  4  [pp.  5-14]  in 
which  a  coherent  receiver  was  assumed,  and  AWGX  used  as  the  interference  model. 
The  three  bits  associated  with  each  symbol  are  identified  as  IB,  2B  and  3B  corresponding 
to  the  most  significant  bit,  the  second-most  significant  bit  and  the  least  significant  bit, 
respectively.  The  results  of  the  BER  analysis  showed  that  the  individual  bit  error  prob- 
abilities were  not  equal  for  the  IB,  2B  and  3B.   In  fact  these  probabilities  are  given  by 

PKIB  incorrect)  =  j Q[^^   _^.^    „. .^  ,    ^  .^^   ^.^    ..^  ^^  ,^ 

=  P;'{3B  incorrect} 

and 
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Figure  8.      Signal  Constellation  for  8-PSK  Showing  Variable  Phase  Parameters 
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Figure  9.      Signal  Constellation  for  16-PSK  Showing  Variable  Phase  Parameters 
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Pr{2B  incorrect}  = 

Q 


E 


Ar 


[  COS  a  +  sin  a]    +  Q[ 


A- 


z-  [  cos  a  —  sin  a] 


^Q\      /~v~  ^  ^°^  "^  "*"  ^^^  "'■-'  ^    -\  /~v~  f-  ^°^  "^  ~  ^^^  °''-^ 


The  overall  BER  for  8-PSK  using  the  DBD  receiver  can  be  expressed  as 


(2.2) 


7  1 

p^{e]  =-^Pr{lB  incorrect}  +  -^  Pr{2B  incorrect} 


(2.3) 


For  16-PSK,  there  are  four  bits  associated  with  each  symbol  which  are  identified  as 
IB,  2B,  3B  and  4B  corresponding  to  the  most  significant  bit,  the  second-most  significant 
bit,  the  third-most  significant  bit,  and  the  least  significant  bit,  respectively.  The  error 
probability  associated  with  each  of  these  bits  is  different  as  given  by  Ref  1  [pp.  461-463], 
under  the  assumption  of  a  coherent  receiver  and  an  AWGX  interference  model.  The 
bit  error  probabilities  for  the  individual  bits  in  16-PSK  are  given  by 


Pr{\B  incorrect}  = 

Q 


4 


2£ 


.V   A( 


cos  a    +  (? 


-i 


Q 


2£ 


a; 


sin  a    +  Q[ 


=  Pr{4B  incorrect} 


2E 


cos  /? 


2£ 


sin^ 


(2.4) 


Pr{2B  incorrect}  =  -z-  Q 


E 


A 


T-  [  cos  a  —  sin  a]    +  Q\ 


t-e( 


a; 


[  cos  a  —  sin  a] 


a; 


[cos/?-sin/?]    +e 


a; 


[  cos  /?  -  sin  /?] 


_E 
A 


r-  [  COS  a  +  sin  a]    • 


E 


n.5) 


A 


e  [  cos  ^  +  sin  J?]    . 
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Pr{3B  incorrect}  = 

4 

2 


+ 


(=5 


/_j\/  ttA'o  Jq 


exp(^  - -^/^(r,  a,  ^)j-(r,  a,  ^) . 


-K\/¥^'^^'"'^^)k 


exp(^  -  ^y;'(r,  a,  PYjgiiT,  a,  j?) . 


t-^ 


2£ 


Ar 


^,(r,  a,  i5)  \dr 


(2.6) 


where  the  integral  in  Eq.  2.6  has  not  been  solved  analytically,  but  has  been  evaluated 
numerically  in  Ref.  1  [pp.  466-467],  and  the  functions /(F,  a,  /?)  and  g,{r,  a,  p)  are  de- 
fined in  Table  1 . 


Table  1. 

FUNCTION  DEFINITIONS 

Index 
i 

/(r,  a,  ^) 

^,(r,  a.  /?) 

1 

sin(r  —  a) 

cos(r  —  a) 

2 

cos(r  -  a) 

sinfr  —  a) 

J) 

cos(r  +  a) 

sin(r  +  a) 

4 

s'mir  +  a.) 

cos(r  +  a) 

5 

simr-/?) 

cos(r-/)) 

6 

cos(r-/?) 

sin(r  -  /?) 

7 

cos(r  +  /?) 

sinir  +  /?) 

8 

sin(r  +  /?) 

cos(r  +  /?) 

The  overall  BER  for  16-PSK  can  then  be  expressed  as 


-r-  ?r{lB  incorrect}  -t-  —  P?'{2B  incorrect}  +  —  Pr{3B  incorrect} 


(2.7) 
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These  results  have  been  presented  here  so  that  results  in  the  sequel  can  be  compared  to 
the  BER  results  associated  with  DBD  coherent  receivers  for  8-PSK  and  16-PSK  signal- 


ing. 
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III.     MFSK/QPSK 

A.     MFSK/QPSK  AS  A  BIORTHOGONAL  SIGNAL  SET 

If  the  features  of  M-ar\-  FSK  and  QPSK  modulation  are  combined,  a  hybrid  digital 
modulation  scheme  results  where  signals  are  both  phase  modulated  and  frequency 
shifted.  This  new  modulation  technique  will  have  four  signals  at  each  of  the  operating 
frequencies  corresponding  to  the  four  different  possible  phase  modulations.  The  signal 
set  can  be  expressed  as 


s,{t)  =  ^J^  cos[2nf,i  +  if:\ 


(3.1) 


for  i  =  0,  1,  2,  3   and  0  <  r  <  r,  ,  and 


s,{t)=J^  cos[27r/2/  +  (/-4)-^] 


(3.2) 


for  i  =  4,  5,  6,  7   and  0  <  r  <  T, . 

Integrating  the  square  of  the  signal  So{[)  over  the  interval  0  <  [  <T,  yields 


rz 


rr^ 


SQ{t)dt  = 


^^cos\2nf,iyt 


=  £ 


1  + 


sin[47r/;rj 


(3.3) 


If  the  product  f^T,  is  an  integer,  or  iffiT,  >  >  1  ,  then 


rr. 


sl{t)dt  =  E 


(3.4) 


which  corresponds  to  the  signal  energy.    Similar  integration  yields  identical  results  for 
5,(0  ,i=l,  2,  ...,  7. 
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Applying  the  Gram-Schmidt  orthonormalization  procedure  Ref.  5  [pp.  266-269]  to 
the  signal  set  described  by  Eqs.  3.1  and  3.2  generates  a  set  of  orthonormal  basis  func- 
tions {(/)//)}  having  the  property 


^■^  \U     for  i=j\ 

0;     for  li-j 


(3.5) 


so  that  the  signals  s,{t)  can  be  synthesized  using  a  linear  combination  of  basis  functions, 


I.e., 


N 


40  =  Yj^ij<f>p) ;    i  =  o,i,  ...,7 


y=i 


where  the  5„  coefficients  are  defmed  as 


(3.6) 


5y     = 


z 


s-^t)(i)h)dt 


(3.7) 


The  first  basis  function  is  obtained  by  normalizing  5o(/)  .   From  Eq.  3.4,  we  have 


</>iW 


^o(0 


\ 


'£ 


—  cos[2;:/iO 


(3.8) 


so  that  5o(/)  =  ^ t  (f)i{t),  and  the  coefficient  %  =  JE  .    From  Eq.  3.7,  5ii  is  given  by 


^11  = 


S^{l)(j)^{[)dt 


^J^  cosi2nft  +  f)J^  cos{2nft)dt 


^\  \  cos[  f  ]  +  cos[2rr(2/,)r  +  f  ]}.//  =  0 


(3.9) 


Defining  the  auxiliary'  function  yi{t)  as 


yiW  =  5,(/)-5,i0i(/) 


(3.10) 
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where  the  second  equality  follows  from  Eq.  3.9,  the  second  basis  function  can  be  ob- 
tained by  normalizing  the  auxiliary  function.   That  is 


HO  = 


yi(0       s^U) 


E 


—  cos(27r/,/  + y) 


Using  next  the  coefficients  53,  and  522  to  obtain  the  auxiliary'  function   y2(/)  as 


where 


and 


V2U)  =  -52(0  -S2i(f>i{t)-  S224>i{t) 


h\ 


=      's2{i)(i)^{i)dt 

/^  cos(2;r/;/  +  7r)     /—  cosi2nf^t)dt 

\  cos[7r]  +  cos[2;r(2/;)f  +  ny^di  =  -  Je 


^      •'0 


(3.11) 


(3.12) 


(3.13) 


522  = 


r. 


S2{t)4>lWt 


0 


!  I I 

^/^  cos(27r/;r+7r)^-^  cos{2nf,t  +  ^)dt 


'E 
T 


rr. 


37Z 


{  cos[  y  ]  +  cos[27:(2/i)r  +  -^  ]} J/  =  0 


(3.14) 


results  in  Eq.  3.12  becoming  equal  to 

y2(/)  =  52(/)  +  7£(/)i(/) 


y-  cos(27r/;r  +  7r)  +  ^  /^  cos(27r/;/) 


2£ 


r. 


[-cos(27r/;0  +  cos(2;r/;/)]  =  0 


(3.15) 


This  indicates  that  52(0  is  completely  described  as 
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S2{t)  =  -jE(i)^{t) 


(3.16) 


and  thus  no  new  basis  function  is  introduced  by  S2{i) 
The  coefficients  for  s^it)  are  given  by 


^31 


-^  cos{2nf,t  +  ^)J^  cos{2nf,t)dt 


iL 

T. 


L 


{  cos[  -^  ]  +  cos[2<2/;)/  +  -^  ^]di  =  0 


(3.17) 


and 


0 


'2E 


jrr 


n 


cos{2nf^i  +  -^)^/-j:-  cos{2nf^t  +  —  )dt 


IE 


Ts 


T, 


'{  cos[7r]  +  cos[2;r(2/,)0)^^  =  -  JE 


(3.18) 


so  that  the  auxiliary  function  y^{i)  becomes 


V'sCO  =  ^3(0  -  ■S32</>2(0 


n 


^  cos{2nf,:+3f)  +  J^  cos{2nf,t  +  f) 


^  [  -  cosi2nf,t  +  Y )  +  cos{2nf,t  +  y  )]  =  0 


(3.19) 


Thus,  as  with  52(0,  s^it)  introduces  no  new  basis  functions.    Continuing  this  procedure 
for  the  signal  s^{t)  ,  we  have 
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Too 

■S41  =  54(00  i(/M/ 

*'-oo 


2£  /  2 

—  cos{27:f2t)/—  cos{2nf^t)d[  (3.2O) 


2V£    fT-. 


cosilTzfjt)  cos{2nf^i)cli  =  0 


where  the  last  equality  is  due  to  the  assumption  that  the  frequencies/  and  j^  are  spread 
far  enough  apart  to  maintain  orthogonality  as  in  binary  FSK  modulation.   Similarly 


;42  =       s^{l)4)2{t)dt 


VT"  ^°^^^''f2t)-J—  cos{2nf^t  +  —  )dt  (3.21) 

'cos(2;r/2/)  cos{2nf^t  +  -y  ^r  =  0 


■'o 

27^"    C^s 


T 


so  that  the  auxiliarv*  function  yj,i)  becomes 


y^t)  =  s^{t)  -  Yj^Mt)  =  s^{t)  (3.22) 


7=1 
Normalizing  yjj)  in  order  to  obtain  the  third  basis  function,  results  in 


Ut)  =  ^  =  .  /"F  cos(2;r/2/)  (3.23) 


Next,  evaluating  55^  for  j  =  1,  2,  3   yields 


•551=      's^{i)4>^{t)dt 


—  cos(27r/2r  +  y)     /—  cos{2nf^t)dt  (3.24) 

■'0 

27^  r^ 


T 

^s      •'0 


I    'cos(27r/;r  +  ^)cos(27r/ir)t//  =  0 
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and 


•S52  =      'ssir)({)2{t)dt 
•'0 

J'Y'  co^(27r/2r  +  Y)    /—  cos(27r/;r  +  yyr  (3.25) 


''0 

2  ^-■^^^'"^'^1'       2 


=  -^^  [  ''cos(27:/2r  +  -f- )  cos(27r/;r  +  -^  ^r  =  0 


where  again,  the  frequencies/  and  j^  are  chosen  to  maintain  orthogonality  as  explained 
above.    Furthermore, 


■^53  =      '55(003(0^^ 


'T, 


-y-  cos(27r/2/  +  Y)     /y-  cos(2nf2t)dt  (3.26) 

■'0 

''  "^^i  cos[  ^  ]  +  cos[27r(2/2)/  +  f  ]}^/  =  0 


Ts 


0 


so  that  with  all  three  of  the  coefficients  equal  to  zero,  5j(f)  must  be  orthogonal  to  all  of 
the  basis  functions  generated  thus  far.  The  fourth  basis  function  is  then  obtained  by 
normalizing  s,{t)  ,  that  is 


Ut)  =  ^=  ^/-f  cos(2^/2^  +  f )  (3-27) 


\ 


It  can  be  easily  shown  that  the  remaining  signals  S(,{i)  and  s-:{i)  can  be  expressed  as 

57(0  =  -  JE  (P^iO 

Thus  the  coefficients  5^3  and  5.4  are  both  equal  to  —  JE  .  Since  the  set  {0//)}  represents 
an  orthonormal  basis  in  a  vector  space,  the  signals  s,{t)  can  now  be  interpreted  geomet- 
rically as  vectors  having  components  equal  to  the  5,^,  namely 

5,(0    -»•   5  =  (5,1,  5,2,  . . .  ,  5,v) ;     /  =  0,  1,  .  . .  ,  A/  -  1  (3.29) 
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where  X  is  the  number  of  linearly  independent  orthormal  basis  vectors,  and  M  is  the 
number  of  signal  vectors. l  It  is  well  known  that 


.V 


sf{t)dt=  )  (5y)'=  |^.|^  =  £,-  (3.30) 


»/ 


0 


;=i 


where  £,  is  the  energy  of  the  ith  signal.    For  the  case  of  2FSK/QPSK,  M  is  equal  to  8, 
N  is  equal  to  4,  and  the  vectors  representing  the  signals  are 

?o  =  (  +  v^,_0,  0,  0) 

?,  =  (0,  +  vr,  0, 0) 

1^^{-Je,o,  0,  0) 

54 = (0, 0,  +  vr ,  0) 

Ts  =  (0,  0,  0,  +  ^lE) 
7,  =  (0,  0,  -  v£  ,_0) 
7-  =  (0,  0,  0,    -  V£  ) 

Since  there  are  eight  different  possible  transmitted  signals,  each  signal  can  represent 
three  bits  of  information. 

The  average  energy  and  also  the  dimensionahty  of  the  signal  set  can  be  minimized, 
bv  subtracting  from  each  signal  a  constant  vector  a  such  that 


M-l 

Y^P{7^]\7^-a\'  (3.32) 

is  minimized,  where  P{s,]  is  the  a  priori  probability  than  the  signal  s^t)  has  been  trans- 
mitted.  Such  vector  a  can  be  obtained  from 


M-] 

a 


ivi  —  1 

=  Y,P{7-;}7-,^E{7}  (3.33) 


i=0 


1  The  general  practice  of  labeling  the  number  of  signals  as  M  is  followed  here,  but  the  letter 
M  in  "MFSK/QPSK"  represents  the  number  of  frequency  shifts. 
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For  a  signal  set  of  the  type  described  above,  assuming  all  a  priori  probabilities  P{'s}  are 
equal,  Eq.  3.33  for  a  reduces  to 

a=\Yji  =  ^  (3.34) 

so  that  there  is  no  minimum  energy  transformation  that  would  reduce  the  dimensionality 
or  the  average  energy  of  the  signal  set. 

Figure  10  on  page  28  shows  the  vector  representation  of  the  2FSK/QPSK  signal 
set  shown  as  eight  vectors  in  a  four  dimensional  vector  space.  A  signal  set  consisting 
of  N  orthogonal  signals  augmented  with  the  negative  of  each  signal  is  referred  to  as  a 
biorthogonal  signal  set.  As  shown  by  Figure  10,  the  2FSK/QPSK  signals  form  this  type 
of  set.  For  a  general  biorthogonal  signal  set,  the  number  of  dimensions  of  the  vector 
space  is  always  equal  to  half  the  number  of  signal  vectors.   That  is 

M=2N  (3.35) 

while  the  number  k  of  bits  represented  by  each  signal  is  given  by 

k  =  Xog^M  (3.36) 

Finally,  the  number  Kj  of  signal  frequencies  in  the  MFSK'QPSK  scheme  is  one-fourth 
the  number  of  signals  or 

M=AKf  (3.37) 

The  best  assignment  of  bits  to  the  M  signals  of  the  biorthogonal  signal  set  is  to 
match  those  signals  separated  by  the  greatest  Euclidean  distance  in  the  vector  space, 
with  symbols  having  the  greatest  Hamming  distance.  With  this  rule  in  mind  the  bit  as- 
signments become  obvious.  That  is,  signals  which  are  antipodal  are  assigned  symbols 
which  correspond  to  binary  complements.  Once  this  assignment  is  done,  all  other  vec- 
tors' Euclidean  distances  are  equal,  so  that  any  other  assignment  patterns  can  be  freely 
chosen.   A  typical  bit  assignment  is  shown  for  the  case  of  A/=  8  in  Figure  10. 

B.     OPTIMUM  RECEIVER 

With  the  basis  functions  having  been  specified,  available  analytical  results  on  the 
structure  and  performance  of  receivers  for  discriminating  biorthogonal  signals  received 
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Figure  10.       Vector  Representation  of  2FSK/QPSK 
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in  additive  white  Gaussian  noise  make  it  possible  to  design  an  optimum  receiver  for  the 
discrimination  of  MFSKQPSK  signals.  The  received  signal  r(t)  is  assumed  to  consist 
of  a  signal  component  5,(0  P^^s  a  sample  function  of  a  white  Gaussian  noise  process  n(t). 
That  is, 

r(r)  =  5,(r)  +  nW  (3.38) 

As  shown  by  Ref  5  [pp.  211-235]  the  optimum  receiver  takes  on  the  structure  of  a  cor- 
relation receiver.  A  block  diagram  of  the  general  correlation  receiver  is  shown  in 
Figure  1 1  on  page  30.  The  signal  r(t)  is  processed  by  a  bank  of  correlators,  where  the 
number  of  correlators  in  the  receiver  is  equal  to  N,  the  number  of  dimensions  in  the 
signal  vector  space  model.  The  effect  of  the  additive  white  Gaussian  noise  is  to  add  to 
the  signal  vector  a  noise  vector  ii  having  X  independent  components,  each  of  which  is 
a  zero  mean,  Gaussian  random  variable  with  variance  —r-  .  Thus  the  correlator  outputs 
Tj  form  the  components  of  a  random  vector  r,  i.e., 

7  =  5;  +  n  =  (ri,  r2,  ...,  r^)  (3.39) 

After  generating  the  vector  f  ,  the  components  r^  are  processed  through  a  weighting 
matrix  to  generate  the  dot  products  d,  which  are  given  by 

N 

^i  =  'r  *^i  =  J]rjSy  (3.40) 

The  d,  are  then  summed  with  bias  terms  c,  where 

c,  =  ^(Aoln[P{m,}]-  i:^!');      /  =  0.  1.  ...,  M-l  (3.41) 

Finally,  all  summer  outputs  1,  are  routed  to  circuitry'  which  decides  which  signal  was  sent 
according  to  the  rule: 

'      *  (3.42) 

Then  decide  s-^t)  was  the  transmitted  signal. 

For  the  signal  vectors  of  Eq.  3.31,  assuming  P{s,{t)]  =-77  for  all  /,  the  structure  of 
the  optimum  receiver  simplifies  from  that  of  the  general  correlation  receiver  to  the 
2FSK  QPSK  receiver  shown  in  Figure  12  on  page  32.  The  2FSK/QPSK  scheme  can  be 
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Figure  11.       Diagram  of  a  General  Correlation  Receiver 
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extended  beyond  the  case  considered  by  increasing  the  number  of  signal  frequencies. 
For  example,  in  order  to  increase  the  number  of  bits  per  signal  from  three  to  four,  the 
number  of  signals  M  must  be  doubled  to  16.  This  is  accomplished  by  doubling  the 
number  of  frequencies  K^  from  two  to  four,  and  thus  doubling  the  number  of  basis 
functions  from  four  to  eight.  This  change  will  be  reflected  in  the  receiver  structure  by 
a  doubling  of  the  number  of  correlators  to  eight.  The  new  signal  set  is  also 
biorthogonal,  and  further  extensions  to  MFSK/QPSK  result  in  N  basis  functions,  Kf 
frequencies,  M  signals,  and  k  bits  per  signal.  The  formula  for  the  probability  of  a  correct 
symbol  decision  for  the  biorthogonal  set  is  given  by  Ref  5  [pp.  261-263]  as 


roo  r  Too 

Ps{Q  =  J^  /„(a  -  V £  m  1  -  2j    fMdp 


iV-l 


(3.43) 


where /^(jc)  is  a  Gaussian  distribution  function  expressed  as 

2 


/n(-^)  = 


1 


V 


2no 


exp    - 


2a' 


(3.44) 


and  where 


a  = 


An 


(3.45) 


Using  the  transformation  >•  = 


a 


/£ 


and  rearranging  the  terms  in  Eq.  3.43  yields 
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Using  another  variable  transformation  the  inner  integral  becomes 


(3.46) 


fl 


(3.47) 


Substituting  this  result  into  Eq.  3.46  gives 
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Figure  12.       Diagram  of  Correlation  Receiver  for  2FSK/QPSK 
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/»oo        r  I —       t"^'-' 


PsiQ 


«/ 


l-2Q[y  + 


v£  ^ 


giy)dy  (3.48) 


Thus  we  have  the  probability  of  a  correct  symbol  decision  associated  with  the  optimum 
receiver  expressed  in  terms  of  the  Q  function. 

C.     OPTIMUM  RECEIVER  BER 

In  order  to  compute  the  Bit  Error  Rate  (BER)  of  the  optimum  receiver  shown  in 
Figure  12,  it  is  necessar>"  to  determine  under  what  conditions  a  chosen  bit  will  be  cor- 
rectly indentified  by  the  receiver  for  a  given  transmitted  signal.  The  three  bits  assigned 
to  each  2FSK,'QPSK  signal  are  identified  as  IB,  2B  and  3B  corresponding  to  the  most 
significant  bit,  the  second-most  significant  bit,  and  the  least  significant  bit  respectively. 
We  can  assume  without  any  loss  in  generality  that  5o(r)  is  transmitted.  The  probability 
that  the  IB  is  correctly  identified  at  the  receiver,  conditioned  on  the  assumption  that 
So{t)  was  transmitted  can  be  written  as 

Pr{  1 B  Correct  \  So]  =  Pr{r,  =  max  r,- 1  Sr,]  +  Pr{  — r,  =  max  r.- 1  Sq] 

"  '         y       \_:  ^        J  \^  (3.49) 

+  Pr{r2  =  max  r.- 1  Sq}  +  Pr{  -r^  =  max  r,- 1  Sq] 
j      -'                               j 

Now 

Pr{x^  =max  r.-ITg} 

j  ^     (3.50) 

=  Pr{r^  >  -r,,  r^  >  Vj,  r^  >  -r^,  r,  >  r3,  r,  >  -r^,  r^  >  r^,  r,  >  -y^  \  Sq] 

If  we  further  condition  on  r,,  that  is  r,  =  /?,  ,  since  r^  >  —r^  only  for  i?i  >  0  ,  then  this 
conditioning  will  be 

Pr{r^  —  max  r.lTo,  y^=  R^>  0} 

'■  ^  (3.51) 

=  Pr{R^  >  rj,  R^  >  -Y2,  R\  >  h,  /?,  >  -r3,  R^  >  r^,  R^  >  -r^  I  Sq,  r^  =  R^>  0} 

Observe  that 

R^  >  Tj,  /?,  >  -Tj    =>     -/?!  <  Fy  <  /?,  (3.52) 

so  that 
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Pr{ri  =  max  r.-I^Q,  r^=^  R^>  0} 

^  (3.53) 

=  Pr{  -R^  <  r2  <  /?i,    -/?i  <  r3  <  R^,    -R^  <r^<  R^\7q,  r^  =  R^>  0} 

Assuming  So{[)  was  transmitted,  then  from  Eqs.  3.8  and  3.38  we  have 


r(/)  =  vT(/)i(r)  +  «(0 


(3.54) 


so  that 


=  V^  +  Hi  , 


'1  ~  V 


1  ,     h  =  ^2>    h  =  ^3  >    '■4  =  n4 


(3.55) 


Due  to  the  independence  of  the  random  variables  n^  ,  the  set  {r,}  consists  of  independent 
Gaussian  random  variables  having  probability  density  functions  (p.d.f.'s) 
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2o' 


(3.56) 


and 


/r,(«>) 
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/27ta' 


exp 


-r: 


IG 


;    for;=2,  3,  4 


(3.57) 


So  that  Eq.  3.53  simplifies  to 


Pr{-R^<x-^<R^,    -i?,  <r3<7?i,    -R^<x^<R^\%,  x^  =  R^>^) 


=  ]~[Pr{-/^i<r^.</?,iro,  ri  =  /^i>0} 


7=2 


Furthermore,  from  Eq,  3.57 


(3.58) 
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Pr{  -/?,  <  Xj  <  /?,  I  ?o.  rj  =  /?,  >  0}  = 
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exp       ^     ]dx 


(3.59) 


=  l-22(-^);     for;=2,  3,  4 


Using  this  result  with  Eqs.  3.53  and  3.58  yields 


Pr{ri  =  max  tUq,  r,  =  /?,  >  0}  = 


i^, 


-<^(4^) 


(3.60) 


so  that  integrating  over  all  R^>0  yields 


Pr{ri  =max  rlTg} 
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.--(:■ 


3 


2 


exp 
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1-20J+-V 


^-^''{^y 


(3.61) 


Note  that  this  expression  approaches  unity  as  —j  -^  oo  as  expected.   Consider  now 

Pr{  — r2  =  max  r.-ITo) 

^  (3.62) 

=  Pr{-r2>r^,    -r2>-ri,    -r2  >  rj,    -r2  >  r3,    -rj  > -r3,    -r2  >  r4,    -r2>-r4l5o} 

If  we  further  condition  on  r^,  that  is  Tj  =  /?2  ,  since  — Fj  >  Tj  can  only  be  satisfied  for 
/?2  <  0  ,  then  this  conditioning  will  be 

Pr{  -r2  =  max  r.- 1  Tq,  rj^  R2<  0} 

^'  ^  (3.63) 

=  Pr{r^  <  -R2,  ri  >  R^,  r^  <  -R2,  r^  >  R2,  r^  <  -R2,  r^>  Ril 'sq,  Tj  =  /?2  <  0} 

Let  R'2  =  —  R2  ,  so  that  in  terms  of  R'2  , 
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rj  <  -R2,  rj  >  /?2    =>    Yj  <  R'2,  Tj  >  -R'2 


-R'2<rj<R'2-    y=l,  3,  4(3.64) 


and  the  requirement  /?2  <  0    =>    i^'j  >  0  ,  so 


/•/-{ -r2  =  max  r,|5o,  Y2  =  —R'2,  R'2>0] 


FotJ=  1  , 


(3.65) 


Pr{  -R'2  <  r,  <  ^^'2} 

rR'2 
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(3.67) 


=  i-2e(^) 


so  that 
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Pr{  -r2  =  max  r.ITo,  r2  =  -R'2^  ^'i  >  0} 


^'0 


^2 


-2<4-) 


,_ef£4^Vd 


R'l-.IE 


(3.68) 


Since  Fj  takes  on  values  R2,  we  can  define  a  new  random  variable  r'j  =  — rj  where  r'^  takes 
on  values  R'2,  i.e.,  if  r^  =  R^,  then  r'j  =  —  R2  —  R'2  .   So  the  condition 

r2  =  -i?'2,  i^'2>0   =>    -r'2  =  -i?'2,  /?'2>0   =*   r'2  =  7?'2,  i?'2>0       (3.69) 

and  given  the  symmetric  p.d.f  for  r^  conditioned  on  Jq  ,  we  see  that  the  conditional  p.d.f 
of  r'2  is  no  different  than  the  conditional  p.d.f  of  rj  ,  so 


Pr{  —Xj  =  max  r- 1  Sq]  =  Pr{Y'2  =  max  r- 1 50} 
j  J       ■' 

f*oo_  2r 
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Through  a  transformation  of  variables  this  expression  is  reduced  to 


Pr{  —Y2  =  max  Xj\'sq] 
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(3.71) 

\dy 


Note  that  this  last  expression  approaches  zero  as 


00  .   Consider  now 


Pr{x^  =  max  xA'sq] 

J 

=  Pr{Y^  >  ri,  r3  >  -r,,  r^  >  rj,  r^  >  -r^,  r^  >  -rj,  r^  >  y^,  y^  >  -yJ7q]  (^-^2) 
=  Rr{Y^  <  h^  ri  >  -h^  ^2  <  Y2,  Y2  >  -r3,  2y^  >  0,  r4  <  r3,  y^  >  -y^I^q] 

If  we  condition  on  r^  that  is  y^  =  R^  ,  we  must  have  R^  >  0,  or  we  cannot  satisfy  the 
2r3  >  0  condition.   Thus 

Pr{Y2  =  max  n  |  Sq,  r3  =  i?3  >  0} 

j  (3.73) 

=  Pr{-i?3<r,  <7?3,    -/?3<r2</?3,    -R^  <y^<  R^ITq,  y,  =  Rj>  0] 
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From  the  mathematical  form  of  this  equation,  it  is  obvious  that  this  probability,  when 
averaged  over  all  /?3  >  0  produces  the  same  probability  as  that  given  by  Eq.  3.71,  which 
also  is  identical  to  Pr{  —r^  =  max  rj^o)  .   Thus,  we  obtain  from  Eqs.  3.49,  3.61  and  3.71 

—         3 

\\-2Q{y  +  ^)\g{y)dy 
v^  ^  -■  (3.74) 
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If  the  2B  is  selected  then 


Pr{2B  Correct\'sQ]  =  Pr{r■^  =  max  rylTg)  +  Pr{r2  =  max  r,|5o} 

+  Pr{r2  =  max  r.-  i  Tq)  +  Pr{u  =  rn^x  r,- 1  Sq] 
J  J 

The  first  term  on  the  right  hand  side  of  Eq.  3.75  is  given  by  Eq.  3.61,  while  the  last  three 
terms  are  all  of  the  same  form  as  Prlr^  =  max  rJJo)  from  Eq.  3.72  and  will  all  yield  the 
same  probability  as  given  by  Eq.  3.71.  Therefore,  the  probability  that  the  2B  is  correct 
is  the  identical  to  the  probability  that  the  IB  is  correct.  It  can  be  shown  in  similar 
fashion  that  the  probability  of  correctly  identifying  the  3B  is  also  equal  to  that  for  the 
IB.  Thus  Eq.  3.74  is  also  the  expression  for  Pi,{Correc[\so]  the  average  probability  of 
decoding  a  bit  correctly  given  that  So{t)  was  transmitted. 

Now  consider  the  case  of  a  given  signal  s^{i)  being  transmitted,  where  s^{t)  can  be 
any  one  of  the  signals  sit),  /  =  0,  1,  .  . . ,  7  .  For  any  chosen  bit,  the  probability  that 
this  bit  is  correct  conditioned  on  the  hypothesis  that  s^{[)  was  transmitted  can  be  ex- 
pressed in  a  form 


Pr{Chosen  Bit  Correct  \  5)^  =  Pr{Yi  =  max  r, 1 5);}  +  /^ Prlr^  =  max  r,  I  SjJ     (3.76) 

where  r,  and  the  y„  terms  are  the  signed  correlator  outputs  ±  r,  .  The  output  r,  corre- 
sponds to  5;^  in  that  if  r^  is  the  maximum  of  all  the  ±  Vj  ,  then  the  receiver  decides  sjj) 
was  transmitted.  The  outputs  r„  correspond  to  three  other  s,  which  if  selected  by  the 
receiver,  would  result  in  the  chosen  bit  still  being  correct.    The  first  term  on  the  right 
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hand  side  will  again  yield  the  same  probability  as  that  given  by  Eq.  3.61.  The  terms  in 
the  summation  are  all  of  the  same  mathematical  form  as  Eq.  3.72,  and  so  are  equal  to 
the  probability  given  by  Eq.  3.71.  Thus  Eq.  3.74  is  the  expression  for  the  probability 
that  any  chosen  bit  is  correct  for  any  given  signal  transmitted,  and  therefore  yields  the 
overall  average  probabihty  of  a  correct  bit,  Pi{C]  . 

Some  insight  into  the  total  probability  of  a  correct  bit  decision  Pi{C],  may  be  gained 
by  comparing  Eq.  3.74  with  the  probability  of  a  correct  symbol  decision  Ps{C}  given  by 
Eq.  3.48.  Since  N,  the  number  of  basis  functions  is  equal  to  4,  the  first  term  in  the  ex- 
pression for  PbiC]  is  equal  to  P,{C].  This  should  come  as  no  surprise,  since  the  deriva- 
tion for  the  probabihty  of  the  IB  being  correct  began  with  Eq.  3.43  where  the  first  term 
is  in  fact  the  probababilty  that  this  bit  is  correct  because  the  receiver  made  the  correct 
svmbol  decision. 

The  second  term  in  Eq.  3.74  represents  the  probabihty  that  the  bit  is  correctly 
identified  even  though  the  receiver  has  made  a  symbol  error.  Recall  from  Figure  10  that 
only  those  signals  which  have  vector  representations  orthogonal  to  each  other  have  bits 
in  common,  and  for  2FSK,QPSK  there  are  three  such  signals  meeting  this  requirement 
for  anv  one  fixed  sienal.  The  probabihties  of  erroneouslv  choosine  one  of  these  sienals 
over  say  s^ii)  ,  are  all  equal  as  discussed  above.  Furthermore,  the  conditional  probabihty 
of  (incorrectly)  choosing  any  of  the  signals  orthogonal  to  s^it),  can  be  shown  to  have  the 
mathematical  form  of  Eq.  3.72.  Thus  Eq.  3.71  specifies  the  probability  of  erroneously 
choosing  any  one  of  the  signals  orthogonal  to  5o(r)  . 

The  above  results  for  2FSK,QPSK  can  be  extended  to  MFSK/QPSK,  that  is.  a 
multi-frequency  four-phase  signaling  scheme.  The  structure  of  the  optimum 
MFSKQPSK  receiver  is  the  same  as  that  of  the  general  correlation  receiver  shown  in 
Figure  1 1  on  page  30  where  the  bank  of  A'  correlators  corresponds  to  a  pair  of 
correlators,  one  each  for  the  In-phase  and  Quadrature  channels,  for  each  of  the  K^  fre- 
quencies. The  expression  for  the  conditional  probability  of  a  correct  signal  decision  has 
already  been  given  in  Eq.  3.48.  The  expression  for  the  conditional  probability  of  decid- 
ing in  favor  of  a  signal  that  is  orthogonal  to  the  transmitted  signal,  yet  does  not  result 
in  a  bit  error  takes  on  a  form  similar  to  that  derived  for  the  2FSK/QPSK  case.  To 
demonstrate  this,  assume  that  So{t)  has  been  transmitted  and  examine  the  probability 
that  one  of  the  signals  orthogonal  to  Sq{i)  is  erroneously  choosen  over  So{i)  .  This  corre- 
sponds to  the  probabihty  that  for  some  r^ ,    A  =  2,  3,  .  .  .  ,  iV 
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Pr{r^  =  max  r.-Ug} 

'  (3.77) 

=  Pr{ri,  >  ri,  r^  >  -r^,  r^,  >  Tj,  r^,  >  -Yj,  . . . ,  r^^  >  -r^,,  . . . ,  r^,  >  r^,  r^  >  -r,vUo} 

Further  conditioning  on  r^  yields 


Pr{r^  =  max  Tj  I  Tg,  r^^  =  R;,,  R/,  >  0} 


■S' 


=  Y\M-Rk<rj<Rk\%,  r,  =  R„  R,>()] 


(3.78) 


The  first  term  in  the  summation  (/=  1)  is  the  same  as  that  given  by  Eq.  3.66  for 
2FSK,QPSK.  All  the  other  terms  (/  =  2,  3,  .  .  .  ,  A';  J^  k)  are  given  by  Eq.  3.67,  so 
that 


M^k  =  i^ax  r- 1 5o,  r;,  =  R^,,  Rj,  >  0} 
J 

_,  A-2  r 


=  [-<^) 


/  R'2  +  JE   \         (  R'l-  JE 
1-g        \^        ]-Q        \^ 


(3.79) 


Nov^^  integrating  over  all  r^  >  0  yields 


Pr{rj.  =  max  r.-Uo) 
J 


=     [1 

'o 


J(2C>0] 


A-2 


l-a>-  + 


JE 


n(       ^^ 


J    V 


-=^  exp 
''27r  ^ 


;3.80) 


Since  there  are  M  —  2  sisnals  orthogonal  to  5o(/)  ,  and  for  anv  chosen  bit,  half  of  these 

(M  -  2) 
signals  have  this  chosen  bit  in  common  with  the  one  in  5o(/),  yields =  A  —  1 

ways  of  making  a  symbol  error  while  correctly  indentifying  the  bit  in  question.  Thus  the 

expression  for  the  probability  of  a  correct  bit  decision  conditioned  on  the  signal  sj^t) 

being  transmitted,  is  given  by 
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Pi,{Correct\sQ] 

f*co 


./F 


,V-1 

1^1-2^0-  +  ^)]      g(y)dy 


(3.81) 


j*oo 

+  (,V-1)       ll-2Qiy)-] 


N-2 


(  J^ 


-4-#)' 


giy)dy 


Using  steps  similar  to  those  used  in  the  analysis  of  the  2FSK/QPSK  receiver  as  well  as 
those  presented  above,  it  can  be  shown  that  Eq.  3.81  is  also  the  expression  for  the  overall 
average  probability  of  a  bit  being  correctly  identified  by  the  MFSK/QPSK  optimum  re- 
ceiver Pt,{C]  . 

Since  the  integrals  involved  in  calculating  the  bit  error  probability  for  the 
MFSK  QPSK  optimum  receiver  cannot  easily  be  expressed  in  closed  form,  the  BER  can 
only  be  evaluated  using  numerical  techniques.  Appendix  B  presents  a  computer  pro- 
gram in  FORTRAN  language  used  to  compute  the  actual  bit  error  probability  for  the 
optimum  receiver  by  employing  numerical  integration.  An  upper  bound  for  the  receiv- 
er's bit  error  probability  is  given  by  Ref  6  [pp.  255-257]  as  a  union  bound,  namely 


Pt{e]^Y 


{M-2)0{ 


kE, 


A'r 


+  a 


IkE, 


Ar 


(3.82) 


This  bound  is  a  eood  approximation  to  the  actual  BER  for  A  >  3  at  moderate  to  larse 

A 

D.     DBD  RECEIVER  FOR  MFSK/QPSK 

Another  receiver  structure  for  MFSK/QPSK  signals  is  implementable  by  using  DBD 
techniques.  Ref  2  [pp.  19-25]  A  diagram  of  the  structure  of  this  receiver  is  shown  in 
Figure  13  on  page  42  for  the  2FSK;QPSK  case.  The  DBD  receiver  uses  the  same 
correlator  subsystem  as  the  optimum  receiver  Figure  12  on  page  32,  however  the  two 
receivers  differ  in  the  way  the  correlator  outputs  are  processed.  Whereas  the  optimum 
receiver  processed  these  correlator  outputs  through  a  weighting  network  and  then  se- 
lected its  largest  output,  the  DBD  receiver  routes  the  correlator  outputs  through  a  two 
stage  summing  network  and  then  compares  the  summer  outputs  to  a  zero  threshold  in 
order  to  arrive  at  binarv  decisions. 
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E.     MFSK/QPSK  DBD  RECEIVER  BER 

The  analysis  for  the  derivation  of  the  BER  performance  for  the  2FSK;QPSK  DBD 
receiver  is  similar  to  that  carried  out  for  the  performance  analysis  of  the  optimum 
(symbol  discrimination)  receiver.  We  begin  by  assuming  Sg{[)  is  transmitted.  Since  the 
correlator  outputs  are  identical  for  both  receivers  in  question,  the  analysis  on  the  per- 
formance of  the  DBD  receiver  can  proceed  with  the  correlator  outputs,  which  are  again 
characterized  as  independent  Gaussian  random  variables  with  density  functions  given 
by  Eqs.  3.56  and  3.57.  Note  from  Figure  13  that  there  is  a  path  from  each  correlator 
output  to  the  input  of  all  second  stage  summers.  The  outputs  of  these  summers,  denoted 
b„ ,     «  =  1,  2,  3  ,  are  given  by 

bi  =  ri  -  r2  +  r3  -  r^ 

bj  =  ri  +  r2  +  r3  +  r^  (3.83) 

b3  =  -r,  -  r2  +  r3  +  r4 

Each  b„  will  also  be  a  Gaussian  random  variable  having  a  mean  equal  to  the  sum  of  the 

means  of  the  r,,  and  due  to  the  independence  of  the  r^.  the  variance  of  b„  will  consist  of 

A' 

the  sum  of  the  variances  of  the  r^  ,  each  of  which  equals  -y-.   Thus  we  have 

b,  ~  A-(v  £  ,  2Ao) 

b2  ~  N{JE  ,JNo)  (3-84) 

b3  ~  A-(  -  JE  ,  2Ao) 

where  N{a,  b)  is  used  to  denote  a  Guassian  probability  density  of  mean  a  and  variance 
b.  With  the  statistics  of  the  b„  specified  comparing  these  random  variables  to  a  zero 
threshold  yields  the  following  conditional  probabiUties  of  decoding  the  bits  correctly 
namelv 


Pr{\B  CorrecA'so]  =  Pr{h^  >  0}  = 


/»oo 

1 


Q     V27r(2A^) 


exp 


2(2Ao) 


dx 

(3.85) 


=  21-^/1^     ='-2V2A„ 
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/•oo 


Pr{2B  Correci  \  Tq)  =  Pr{\i2  >  0} 


x/27r(2Ao) 


exp 


2(2Ao) 


dx 


(3.86) 


=  l-e 


2An 


/»0 


Pr{35  CoA-r^crl^o)  =  Pr^  <  0}  = 


1 


'27t(2iYo) 


exp 


(^  +  7^ )' 

2(2Ao) 


dx 


(3.87) 


=  1-0 


2a; 


Since  each  of  the  individual  bits  have  identical  probability  of  correct  decoding,  the  cor- 
responding average  probability  conditioned  on  the  assumption  that  Sq{i)  was  transmitted 
is 


and  consequently  the  conditional  bit  error  probability  is 

P^iErroAl^]  =  1  -  Pt,{Con-ecr\7Q]  =  Q(yJ^] 


(3.88) 


(3.89) 


Consider  next  the  case  in  which  any  one  of  the  s,{t) ,  /  =  0,  1,  . . . ,  7  is  transmitted,  for 
any  of  the  receiver  correlators.  If  the  transmitted  signal  is  matched  to  the  local  signal 
correlator,  its  output  Vj  will  be 


Yj  =  ■^'E  +  lij 


and  for  the  correlator  having  antipodal  local  signal  matching,  the  output  will  be 


(3.90) 


F;  =  -  -^^  +  Uj 


All  other  correlators  will  have  an  output  given  by 


r,  =  n, 


(3.91) 


(3.92) 
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Examining  now  any  one  of  the  second  stage  summer  outputs  b„ ,  if  the  bit  corresponding 
to  b„  for  the  signal  sent  is  a  logical  zero,  then  b„  is  statistically  characterized  by 


b,  ~  MJE  ,  2Ao) 


(3.93) 


In  this  case  the  bit  will  be  correct  if  b„  is  greater  than  the  zero  threshold,  and  the  prob- 
abiUty  of  the  bit  corresponding  to  b„  being  correctly  decoded  is  given  by 


/•oo 


Pr{K  >  0}  = 


V2;r(2Ao) 


exp 


jx  -  JE  y 

2{2No} 


dx 


=  l-e 


E 


V   2A( 


(3.94) 


However,  in  the  case  where  the  bit  is  a  logical  one,  then 


b,  -  iV(  -  ^E  ,  2Ao) 


(3.95) 


Now  the  bit  will  be  correct  if  b„  is  less  than  the  zero  threshold,  and  the  probability  of 
correct  decoding  is  eiven  bv 


1^0 


Pr(b,  <  0}  = 


1 


'2;t(2Ao) 


exp 


jx  +  V  £  y 

2(2Ao) 


dx 


=  1-0 


''    E 


2a; 


(3.96) 


For  both  cases  the  probability  of  a  bit  error  is  the  same,  and  it  is  identical  for  all  b,  re- 
gardless of  which  signal  sXO  is  transmitted.  Thus,  we  can  write  the  total  bit  error  prob- 
ability for  the  2FSK,  QPSK  DBD  receiver  as 


Ptie]  =  Q[ 


2Ao 


(3.97) 


Extending  the  2FSK,'QPSK  results  to  the  general  case  of  MFSK/QPSK  in  which 
there  are  M  signals,  N  receiver  correlators,  k  bits,  and  K^  frequencies  where  the 
biorthogonal  signal  relationships  of  Eqs.  3.35  -  3.37  again  hold,  for  each  of  the  random 
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variables  b^  ,  n=  1,  2,...,  k  ,  the  corresponding  variance  will  be  the  sum  of  the  variances 
of  the  n^,  so 


A' 

VariK)  =  l-f  =  -Y-  =  ^A  (3.98) 

and  again  the  mean  of  the  b„  will  be  ±  Je  .  Following  the  same  steps  as  those  carried 
out  in  obtaining  Eqs.  3.90  -  3.97,  the  total  bit  error  probability  for  MFSK/QPSK  DBD 
receiver  is  eiven  bv 


Pt{e]  =  Q\J-f^  I  (3.99) 

F.     BANDWIDTH  OF  MFSK/QPSK 

The  bandwidth  of  MFSK/QPSK  signals  can  be  derived  from  knowledge  of  QPSK 
and  MFSK  bandwidth  requirements.  Each  carrier  is  switched  at  the  symbol  rate  R, 
where  R^  =  —  .  This  gives  each  carrier  a  two-sided  bandwidth  of  IR,  Hz.  To  maintain 
orthogonality  between  signals  at  different  frequencies  in  MFSK  modulation,  the  fre- 
quencies are  separated  by  R,  Hz.  This  same  relationship  must  also  be  maintained  in 
MFSK 'QPSK.  The  spectrum  of  the  2FSK  QPSK  signal  is  plotted  in  Figure  14  on  page 
47.  From  the  graph,  the  total  null-to-null  bandwidth  of  the  2FSK'QPSK  signal  can  be 
seen  to  be 

W^-hR,  (3.100) 

Proceeding  to  the  family  of  MFSK/QPSK  signals,  Figure  15  on  page  48  shows  the 
spectrum  for  such  signals.  By  examining  this  figure,  the  total  null-to-null  bandwidth  can 
be  shown  to  be 

W={Kf+\)R, 

k-2  (3-101) 

=  (2*  '  +  \)R, 

This  general  result  for  MFSK/QPSK  correctly  yields  the  result  for  the  special  case  of 
2FSK'QPSK  in  which  A  =  3,  namely 

W  =  {2^~^  +\)R,=  ZR,  (3.102) 
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which  agrees  with  Eq.  3.100. 

The  bandwidth  efiiciency  of  a  digital  modulation  scheme  is  given  by  Ref.  6  [pp. 
393-395]  as 

V=-^  (3.103) 

and  is  one  wav  to  evaluate  how  efficientlv  different  modulation  schemes  use  the  available 
channel  bandwidth.  For  M-ar\'  PSK,  bandwidth  efficiency  is  given  by  Ref  3  [pp. 
211-213] 

Vmpsk  "  2R 
kR. 


2R, 


(3.104) 


2 


and  for  M-ar>'  FSK  the  bandwidth  efficiency  is 

R, 


^FSK  — 


(M  +  l)R, 
kR, 

k 


2^+\ 


For  MFSK/QPSK,  using  Eq.  3.101  we  obtain  a  bandwidth  efficiency  given  by 

Rt 

if  ^  +  \)R, 

kR, 

(3.106) 


(2^"^  +  \)R 


2^'~^  +  1 


A  table  of  bandwidth  efficiencies  for  M-ar>-  PSK,  M-ary  FSK,  and  MFSK/QPSK  is 

given  in  Chapter  5. 
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IV.     RAYLEIGH  FADING 

A.  THE  RAYLEIGH  FADING  MODEL 

Often,  when  a  signal  is  transmitted  through  a  channel  in  which  multipath  propa- 
gation and/ or  heav;>'  precipitation  conditions  exist,  fading  occurs  with  a  consequent 
change  in  the  signal  amplitude.  One  way  to  model  this  fading  is  to  redefme  the  signal 
amplitude  A  which  up  until  now  has  been  treated  as  a  constant,  as  a  Rayleigh  distributed 
random  variable  A  with  p.d.f. 

A(^)=-4exp(--^W)  (4.1) 

where  u{-)  is  the  unit  step  function,  and  from  Ref  7  [pp.  110-112],  the  value  of /io  is  re- 
lated to  the  second  moment  of  A  as 

E{X']  =  2AI  (4.2) 

In  order  to  evaluate  the  effect  of  Rayleigh  fading  on  the  BER  for  a  given  modulation 
scheme  and  corresponding  receiver,  the  expression  for  BER  under  nonfading  channel 
conditions  must  be  averaged  over  the  distribution  of  the  signal  amplitude.  This  requires 
expressing  the  bit  error  probability  as  a  function  of  signal  amplitude,  weighting  the 
function  with  the  amplitude  p.d.f  and  integrating  this  product  over  all  values  of  A. 

B.  EFFECT  ON  BER  OF  M-ARY  PSK  USING  DBD  RECEIVER 

For  the  8-PSK  DBD  receiver,  the  (non-fading)  bit  error  probabilities  are  given  by 
Eqs.  2.1  -  2.2.  These  probabilities  must  now  be  expressed  as  functions  of  the  random 
amplitude  A.    From  Eq.  1.1  it  is  apparent  that  the  signal  ampUtude  is  given  by 


A  =  J^  (4.3) 


so  that 


£  =  4^  (4.4) 

Since  the  signal  amplitude  is  now  a  random  variable,  the  signal  energy  also  becomes  a 
random  variable  E,  which  is  expressed  as 
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E  = 


A^r, 


(4.5) 


Substitutine  this  result  into  Eqs.  2.1  through  2.3  yields 

Pr [IB  incorrect}  = 


T^(^y¥^"^)He(^^ 


An 


Sin  a 


(4.6) 


Pr{2B  incorrect}  = 


Q[  A  j-zTT  [  cos  a  +  sin  a]  1  +  ^f  A  /  ^ /.    [  cos  a  —  sin  a] 


2Ar 


2Ao 


—  2Q\  A  /  -  .-    [  cos  a  +  sin  a]  \Q[  A  /  ^        [  cos  a  —  sin  a] 


and  as  before, 


(4.7) 


Pr{3B  incorrect}  =  Pr{lB  incorrect} 


(4.8) 


The  bit  error  probability  under  fading  conditions  for  the  most  significant  bit,  namely  IB, 
is  obtained  by  weighting  the  probability  given  in  Eq.  4.6  by  the  p.d.f  of  the  random 
amphtude  A  and  integrating  over  all  A.   That  is 


^Ray{^^  incorrect}  = 


Pr{\B  mconQci]fjJ^A)dA 


\^[n^  '''')^\'^{N-k ''''') 


fx{A)dA 


(4.9) 


This  integral  can  be  solved  using  the  identity  given  by  Ref.  8  [pp.  12-13]  as 


r 


QWAY^{A)dA  = -^ 


1  + 


2  ^2 
a  ^n 


(4.10) 


Using  Eq.  4.2  and  the  linearity  of  the  expectation  operator  we  can  write  an  expression 
for  the  average  signal  energy  under  fading  conditions,  namely  £(E}  =  E  as 
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F  =  f)  ^ 


Ts  2 

=  tA 


Ts 


(2^o) 


(4.11) 


Using  Eq.  4.10  in  order  to  evaluate  Eq.  4.9  yields  the  bit  error  probability  under  Rayleigh 
fading,  namely 


P/ja^{lB  incorrect}  = 


2 


E  2 

-TT-  COS  a 

•^0 


+  1 


_1_ 

4 


^    sin^a  1     +  I 


A 


(4.12) 


Next,  the  BER  under  fading  conditions  for  the  second-most  significant  bit,  namely  2B, 
is  obtained  by  integrating  the  product  of  the  BER  given  by  Eq.  4.7  with  the  Rayleigh 
p.d.f .  The  contribution  from  the  first  two  terms  of  Eq.  4.7  can  be  derived  using  the 
same  steps  as  those  carried  out  above  for  the  evaluation  of  ^^^^^(IB  incorrect}.  Since  the 
third  term  in  Eq.  4.7  involves  the  product  of  two  Q  functions,  a  new  identity  is  required. 
That  is 


Q{aA)Q{MYAiA)dA  = 


/?tan 


-1 


4 


^ 

a 


2n   //?'  +  4 


a  tan 


a.    -{ : 


Ar 


P 


Ac 


(4.13) 


the  proof  of  which  is  presented  in  Appendix  A.  Using  this  identity  for  the  contribution 
of  the  third  term  of  Eq.  4.7  and  combining  this  with  the  contributions  from  the  first  two 
terms  of  Eq.  4.7  yields 
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P;ja^,(2B  incorrect}  = 

_]_  1 

2  ■ 


(4.14) 


where  u^  and  Wj  are  given  by 


i/i  =  cos  a  +  sin  a 
^2  =  cos  a  —  sin  a 


(4.15) 


The  overall  averase  bit  error  rate  under  fadin?  conditions  can  be  obtained  from 


2  1 

Pb,  Ray{^)  =  —  pRay{'^^  mcorrect}  +  —  PRayi'^^  incorrcct] 


(4.16) 


For  the  16-PSK  DBD  receiver  operating  under  Rayleigh  fading  conditions,  much 
of  the  performance  analysis  can  be  carried  out  using  the  same  procedure  as  the  one  used 
for  the  evaluation  of  the  performance  of  the  8-PSK  DBD  receiver  operating  under  simi- 
lar Rayleigh  fading  conditions.  The  bit  error  probabilities  under  non-fading  conditions 
are  given  by  Eqs.  2.4  -  2.6.  Rewriting  these  equations  in  terms  of  the  random  variable 
A,  for  the  most  significant  bit  we  obtain 


Pr{lB  incorrect)  = 


Q\Kh:^ 


(=1  ■- 


(4.17) 


where  the  variables  v,  are  defined  as 
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V]  =  COS  a 
V2  =  cos  /? 
V3  =  sin  a 
V4  =  sin  p 

and  for  the  second-most  significant  bit 

Pr{2B  incorrect}  = 


Yp{^j4;^^  -  ^(^^/i|^-.)4^^/3 


W2 


where  the  variables  w,  are  given  by 

Wj  =  cos  a  —  sin  a 
W2  =  cos  a  +  sin  a 
VV3  =  cos  ^  —  sin  /? 
W4  =  cos  /?  -f  sin  ^ 

For  the  third-most  significant  bit 
P/-{3B  incorrect)  = 


-     n^^f     (     A^r. 


T-gV^^lo'^-^Af^'^^'  -  ^V^'<^'  ^'  ^^• 


1  J  .      lis 

a; 


^     ^(Aj-^g,(r,  a,  (j))}c{r 


/r  f   «  /  A     /, 


(4.18) 


(4.19) 


(4.20) 


(4.21) 


where  the  fiinctions/(r,  a,  /?)  and  g,(r,  a,  /?)  are  defined  in  Table  1  of  Chapter  2.   For 
the  least  significant  bit,  from  Eq.  2.3  it  is  obvious  that 

P/(4B  incorrect]  =  Pr [IB  incorrect)  (4.22) 
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For  the  most  significant  bit,  the  BER  under  fading  conditions  is  obtained  using  Eq.  4.17 
and  the  identity  given  by  Eq.  4.9,  to  yield 


P;j^^{ IB  incorrect}  =-T- 


8 


(=1 


=  Pnayi'^^  incorrect} 


E       2 


,-1 


+ 1 


(4.23) 


where  the  variables  v,  are  defined  in  Eq.  4.18,  and  the  last  equality  in  Eq.  4.23  follows 
from  Eq.  4.22.  The  BER  under  Rayleigh  fading  conditions  for  the  second-most  signif- 
icant bit  is  obtained  using  the  identities  of  Eqs.  4.9  and  4.13  to  yield 


PRayi^^  incorrect}  =  —  -  "T 


4 


;=1 


^     vvA"'  +  l 


2.Vo  '^' 


w^  tan 


-1 


+ 


w, 


+  2 


X 


W-: 


W2  tan 


-1 


+ 


.i  +  2    -f 
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where  the  w,  terms  are  defined  in  Eq.  4.20.    For  the  third-most  significant  bit,  the  inte- 
grals to  be  evaluated  are  all  of  the  basic  form 
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=  [  '  giV,  a,  p)  - 


3   j2 


X   exp(  —  x'^)dxdr 


(4.25) 


where  £,  is  defined  as 


El=        Pi  + 


2Ar 


(4.26) 


and 


r, 


(4.27) 


The  variable  x  is  defined  as  .r  =  Ac,,  so  that  the  evaluation  of  the  integral  of  the  product 
of  the  first  term  in  the  summation  of  Eq.  4.21  with  the  Rayleigh  p.d.f.  becomes 


A^r 


s    A. 


2a; 


y;tr,  a,  /?)  L?/(r,  a,  p)drf^{A)dA 


_i_ 

8 


_1_ 

4 


r..^ 


2  71  A' 


■^0        A 


dT 


0    J 


(4.28) 


^f 


a; 


^/(r,  a,  P) 


A 


^  y;'(r,  «,  ^)  +  i 


dr 


Multiplying  the  second  term  within  the  summation  of  Eq.  4.21  by  the  Rayleigh  p.d.f.  and 

integrating,  yields 
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where  ^,  is  defined  as 
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Carrying  out  the  integration  with  respect  to  A  first,  while  reahzing  that  u{x  —  At)u{A) 


X 


is  nonzero  onlv  for  A  in  the  range  0  <  A  <  -r-,  vields 
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where  a,  is  defined  as 


«i  =  Pi  + 


2A; 


(4.32) 


so  that  Eq.  4.29  becomes 
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Focusing  on  the  integration  with  respect  to  x  only,  we  have  for  the  first  integral  with 
respect  to  .v, 
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and  for  the  second  integral  with  respect  to  ;c 
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where  y,  is  defined  as 
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Thus,  Eq  4.29  becomes 
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From  our  previous  definitions  we  can  rewrite  the  expressions  for  a,  and  y,  as 
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and 
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Observe  that  the  terms  in  Eq.  4.37  can  be  expressed  in  terms  of  the  ratio  of  average 
signal  energy  to  noise  power  spectral  density  as 
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Using  these  results,  the  first  term  on  the  right  hand  side  of  Eq.  4.37  becomes 
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While  the  second  term  on  the  right  hand  side  of  Eq.  4.37  takes  on  the  form 
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Finally,  from  the  third  term  on  the  right  hand  side  of  Eq.  4.37  takes  on  the  form 
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so  that  the  bit  error  probabUity  for  the  third-most  significant  bit  is  obtained  from  the 
sum  of  Eqs.  4.28  and  4.43  -  4.45  .  These  integrals  have  not  been  expressed  in  closed 
form,  however  corresponding  numerical  evaluations  can  and  have  been  carried  out 
without  difficulty.   The  total  BER  under  Rayleigh  fading  conditions  is  finally  given  by 

Pb,  Ray{e]  =  T  ^Rayi^^  incorrect}  -f  —  PRc,y{2B  incorrect} 

1  (4.46) 

+  -J  PRayi'^B  incorrect} 

Appendix  C  presents  a  listing  of  the  FORTRAN  computer  program  used  to  evaluate  the 
BER  under  Rayleigh  fading  conditions  for  the  16-PSK  DBD  receiver.  Data  from  this 
program  is  plotted  and  presented  in  Chapter  5. 

C.     EFFECT  ON  BER  OF  MFSK/QPSK  USING  DBD  RECEIVER 

Recall  that  for  the  MFSK/QPSK  DBD  receiver,  the  BER  is  given  by  Eq.  3.99  as 

^'{^'  =  e(^)  (4.47) 

This  expression  must  be  rewritten  as  a  function  of  the  amplitude  A  in  order  to  determine 
the  effect  of  Rayleigh  fading  on  receiver  performance.  Substituting  Eq.  4.4  into  Eq.  4.47 
yields 


PM  =  Q 

(4.48) 


Weighting  this  probability  by  the  Rayleigh  p.d.f  and  integrating  over  the  variable  A,  the 
average  bit  error  probability  in  the  presence  of  signal  fading  is 
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This  integral  can  be  evaluated  using  the  identity  of  Eq.  4.9  as 
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D.     EFFECT  ON  BER  OF  MFSK/QPSK  USING  THE  OPTIMUM  RECEIVER 

Evaluation  of  the  (non-fading)  BER  for  MFSK/QPSK  using  the  optimum  receiver 
has  been  carried  out  in  Chapter  3,  Section  C.  Yet  a  different  approach  to  determining 
the  BER  under  Rayeigh  fading  conditions  becomes  more  practical  than  using  modifica- 
tions to  the  approach  previously  used  to  evaluate  the  performance  of  the  DBD  receiver. 
Since  the  BER  given  in  Eq.  3.81  as  a  function  of  5,  the  ratio  of  signal  energj^  to  noise 
variance,  namely 
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where 
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has  been  previously  evaluated,  a  transformation  of  variables  is  suggested  to  yield 
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where  the  range  of  s  is  0  <  s  <  oo.   Since  the  transformation  is  linear,  the  p.d.f.  of  s  can 
be  written  as 
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where  SI  is  defined  by 
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The  BER  under  fading  conditions  can  be  evaluated  by  averaging  the  BER  given  by  Eq. 
3. 81  over  the  range  of  the  new  random  variable  s  as 
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\u{s)dycls 


This  integration  also  requires  numerical  techniques,  however  the  program  presented  in 
Appendix  B  can  be  used  as  a  subroutine  for  the  program  listed  in  Appendix  D  which  is 
used  to  evaluate  Eq.  4.56.  Data  from  the  program  of  Appendix  D  is  plotted  and  pre- 
sented in  Chapter  5. 
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V.     RESULTS 

A.     PERFORMANCE  COMPARISONS 

The  symbol  error  rate  (SER)  of  the  receiver  for  M-ary  PSK  signals  is  given  by  Eq. 
1.4,  while  Eq.  1.8  is  used  to  evaluate  the  SER  of  the  receiver  for  M-ary  FSK  signals. 
The  SER  for  MFSK,  QPSK  signals  when  using  the  optimum  receiver  of  Figure  12  on 
page  32  is  given  bv  Eq.  3.48.  Table  2  lists  the  ratio  of  bit  energv  to  noise  power  spectral 
density  (PSD)  level,  -r:-,  required  to  achieve  symbol  error  rates  of  10"^  and  10"^  for  M- 
ary  PSK,  M-ar>-  FSK,  and  MFSK,  QPSK  signaling  schemes  with  k  equal  to  3  and  4 
bits/symbol.    The  table  shows  the  MFSK  QPSK  signaling  scheme  yields  slightly  better 

£ 

performance  than  M-ar\'  FSK  with  a  nominal  0.2  dB  difference,  in  terms  of-rr-  ,  while 

E  ° 

outperforming  M-ary  PSK  by  an  average  -tt-  difference  of  5.2  dB  for  /:  =  3  and  10.4  dB 


A^ 


for  k  =  4. 


Table  2.     PERFORMANCE    COMPARISONS     AT    KEY    SYMBOL    ERROR 
RATES 


Modulation 

Bits  per 
Symbol 

-TT-  (dB)  Required  to  Achieve 
10-3  SER            10-^  SER 

8-PSK 

3 

10.8 

14.4 

16-PSK 

4 

15.5 

1S.4 

8-FSK 

3 

06.5 

08.6 

16-FSK 

4 

05.6 

07.8 

2FSK  QPSK 

^ 
J 

06.2 

08.5 

4FSK  QPSK 

4 

05.5 

07.6 

The  BER  for  MFSK, QPSK  signaling  using  the  optimum  receiver  shown  in 
Figure  11  on  page  30  is  evaluated  using  Eq.  3.81.  Appendix  B  presents  the  FORTRAN 
computer  program  used  to  compute  this  optimum  receiver's  bit  error  probability  by 
employing  numerical  integration  in  order  to  evaluate  Eq.  3.81.  That  is,  Pb{e]  along  with 
the  union  bound  approximation  to  this  probability  as  given  by  Eq.  3.82  is  computed  for 
varying  values  of-r^.   The  program  is  also  used  to  evaluate  Pi,{e]  for  the  .VI FSK  QPSK 
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DBD  receiver,  as  given  by  Eq.  3.99.  Data  generated  by  this  program  is  plotted  in 
Figure  \6  [or  k  =  3,  Figure  17  for  A:  =  4,  Figure  18  for  A'  =  5,  and  Figure  19  for  A  =  6  . 
In  all  of  the  above  referenced  figures,  numerical  values  for  the  BER  obtained  using 
the  union  bound  approximation  yield  plots  almost  coinciding  with  the  actual  BER  plots 
for  the  optimum  receiver  at  moderate  to  high  -rp  values.  This  fact  helps  to  confirm  the 
appropriateness  of  the  results  of  the  numerical  integration  used  to  evaluate  the  actual 
receiver  BER.  From  these  plots,  it  becomes  evident  that  the  xMFSK/QPSK  DBD  re- 
ceiver suffers  a  performance  penalty  when  compared  to  the  performance  of  the  optimum 
receiver.  Table  3  lists  -rr-  ratios  required  to  achieve  BER's  of  10"^  and  10~*  for  both  the 
optimum  and  the  DBD  receiver  structures  for  MFSK/QPSK  modulation.  The  optimum 
receiver  has  an  average  2.5  dB  better  performance  than  the  DBD  receiver  for  /c  =  3 
bits/symbol,  which  increases  to  an  average  value  of  5  dB  for  A:  =  4  in  terms  of-rr-  dif- 
ference.  For  k  =  6,  the  optimum  receiver  has  an  average  10  dB  advantage  over  the  DBD 
receiver. 


Table  3.     PERFORMANCE  COMPARISONS  OF  MFSK/QPSK  RECEIVERS 


Modulation; 
Receivers 

Bits  per 
Symbol 

El. 
-r^  (dB)  Required  to  Achieve 

10-^  BER            10-^  BER 

2FSK  QPSK 

-^ 

:) 

Opt.  Rcvr 

05.8 

09.1 

DBD  Rcvr 

os.o 

11. S 

4FSK  QPSK 

4 

Opt.  Rcvr 

05.1 

08.1 

DBD  Rcvr 

09.8 

13.6 

8FSK  QPSK 

5 

Opt.  Rcvr 

04.5 

07.4 

DBD  Rcvr 

11.8 

15.6 

8FSK  QPSK 

6 

Opt.  Rcvr 

04.1 

06.8 

DBD  Rcvr 

14.0 

17.8 

The  reasons  for  this  increased  performance  penalty  can  be  explained  as  follows. 

Ef, 
First,  as  can  be  seen  by  examining  the  DBD  receiver  BER  curves  as  a  function  of-r^. 

the  performance  of  the  DBD  receiver  degrades  as  k,  the  number  of  bits  per  symbol  is 
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EB/NO  (DD) 
Figure  16.      2FSK/QPSK  Bit  Error  Performance  (k  =  3  bits/symbol) 
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EB/NO  (DB) 

Figure  17.      4FSK/QPSK  Bit  Error  Performance  (k  =  4  bits/symbol) 


67 


5  10 

EB/NO  (DB) 

Figure  18.      8FSK/QPSK  Bit  Error  Performance  (k  =  5  bits/symbol) 
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Figure  19. 


10 

EB/NO  (DB) 

16FSK/QPSK  Bit  Error  Performance  (k 


6  bits/symbol) 
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increased.  Recall  from  the  BER  analysis  for  the  DBD  receiver,  that  the  variance  for  the 
summation  output  b„  is  the  sum  of  the  variances  of  the  correlator  outputs  r^.  When  the 
number  of  bits  (  and  consequenctly  the  number  of  correlators  )  is  increased,  there  is  a 
linear  increase  in  the  variance  of  b„  .  This  variance  increase,  v^'hich  is  strictly  noise  de- 
pendent, causes  the  bit  error  probabihty  increase  also.  On  the  other  hand,  examination 
of  the  optimum  receiver  bit  error  curves  reveals  that,  with  increasing  k,  the  optimum 
receiver  exhibits  improved  performance.  This  improvement  reflects  the  tradeoff  that 
exists  between  bit  error  rate  performance  and  signal  bandwidth,  as  increasing  k  results 
in  a  greater  signal  bandwidth.  Note  however  that  this  improved  receiver  performance 
would  not  be  apparent  if  the  bit  error  probabihty  were  plotted  against  unnormalized 
SNR  instead  of  bit-normalized  SNR,  that  is,  -rp  Ref.  6  [pp.  172-175]. 

B.     BANDWIDTH  EFFICIENCY 

The  bandwidth  efficiency  for  M-arv'  PSK,  M-ary  FSK,  and  MFSK/QPSK  signahng 
is  given  by  Eqs.  3.105,  3.106  and  3.107,  respectively.  Bandwidth  efficiencies  given  in 
bits'Hz  for  values  of  k  ranging  from  2  to  6  bits/symbol  are  listed  in  Table  4  for  these 
three  signaling  schemes.  The  table  shows  that  M-ar\'  PSK  is  more  bandwidth  effiicient 
than  both  the  M-FSK  and  MFSK  QPSK  schemes.  MFSK  QPSK  is  from  three  to  four 
times  as  effiicient  as  M-ar\-  FSK,  with  the  greater  differences  occuring  at  higher  values 
of  k.  MFSK  QPSK  operates  at  67%  the  efficiency  of  M-ary  PSK  for  A  =  3  bits/symbol, 
and  at  about  40^/ o  the  efficiency  of  M-ar\-  PSK  for  /c  =  4  bits  svmbol. 


Table  4.     BANDWIDTH  EFFICIENCIES  FOR  MARY  PSK,  MARY  FSK,  AND 
MFSK/QPSK 


Number 
of  Bits 

Number 
of  Signals 

Bandwidth  Efficiency 
MPSK                  MFSK             MFSK  QPSK 

4 

1.00 

0.40 

— 

3 

8 

1.50 

0.33 

1.00 

4 

16 

2.00 

0.24 

0.80 

5 

32 

2.50 

0.15 

0.56 

6 

64 

3.00 

0.09 

0.35 

C.     EFFECTS  OF  RAYLEIGH  FADING 

The  BER  of  the  8-PSK  DBD  receiver  operating  under  Rayleigh  fading  signal  con- 

Eu 


ditions  is  given  by  Eq.  4.16  as  a  function  of  average 


A'r 


This  equation  has  been  used 
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£, 


to  evaluate  receiver  performance  as  a  function  of  average  -rr-  for  various  values  of  the 
angle  alpha.  The  plots  shown  in  Figure  20  on  page  72  demonstrate  that  the  best  re- 
ceiver performance  occurs  for  a  =  22.5°.  The  curves  for  a  =  20.0°  and  a  =  25.0°  and  for 
a  =  17,5'  and  a.  =  27.5°  are  coincident,  and  therefore  appear  as  one  line  for  each  pair. 

Appendix  C  presents  the  FORTRAN  program  used  to  evaluate  via  numerical  inte- 
gration the  bit  error  probability  under  Rayleigh  fading  signal  conditions  for  the  16-PSK 
DBD  receiver.  Data  generated  by  this  program  is  plotted  in  Figure  21  on  page  73  for 
various  values  of  the  angle  a  while,  maintaining  the  relationship  /J  =  45°  —  a.  This  plot 
shows  the  best  receiver  performance  is  achieved  for  a  =  11.25°  .  Curves  for  a  =  9.75° 
and  a  =  12.75°  as  well  as  for  a  =  8.25°  and  a.  =  14.25°  are  coincident  and  appear  as  only 
one  line. 

The  receiver's  BER  for  MFSK  QPSK  using  the  DBD  scheme,  is  given  by  Eq.  4.50 
when  Rayleigh  fading  signal  conditions  exist.  Figure  22  on  page  74  illustrates  the  plot 
of  this  BER  for  values  of  k  ranging  from  3  to  5  bits/symbol.  The  best  performance  is 
obtained  for  k  equal  to  3,  and  performance  degrades  with  increasing  values  of  k.  The 
performance  curve  has  a  near  linear  characteristic  for  values  of  average  -—  exceeding 
about  5  dB  as  opposed  to  the  "waterfall"  characteristic  of  the  BER  curves  when  Rayleigh 
fading  conditions  were  not  present. 

The  BER  for  MFSK  QPSK  signaling  using  the  optimum  receiver  under  Rayleigh 
fading  signal  conditions  is  given  by  Eq.  4.56  .  Appendix  D  presents  the  listing  of  the 
FORTRAN  program  used  to  evaluate  this  BER,  which  involves  double  numerical  inte- 
gration. Data  generated  by  this  program  is  plotted  in  Figure  23  on  page  75  for  values 
of  k  ranging  from  3  to  5  bits 'symbol.  The  curves  show  a  somewhat  surprising  result  in 
that  the  best  performance  is  obtained  when  k  is  equal  to  3,  and  performance  degrades 
with  increasing  values  of  k.  This  result  is  opposite  to  that  encountered  when  Rayleigh 
fading  conditions  were  not  present,  in  which  increasing  k  produced  improved  receiver 
performance.  The  curves  of  Figure  23  also  exhibit  nearly  linear  characteristics  for  \  al- 
ues  of  average  -rr-  greater  than  approximately  0  dB. 

Table  5  on  page  76  presents  a  summary  of  the  average  -r^  values  required  to 
achieve  bit  error  rates  of  10"'  and  10-^  for  all  of  the  receivers  discussed  above,  under  both 
non-fading  and  Rayleigh  fading  signal  conditions.  For  the  MFSK /QPSK  receivers,  in 
order  to  simplifiy  the  comparisons,  only  the  results  for  k  =  3  are  presented.  For  the 
MFSK  QPSK  receivers,  the  performance  degradation  is  about  the  same  for  both  the 
optimum  and  the  DBD  receivers.  In  order  to  achieve  a  BER  of  10"'  both  receivers  suffer 
a  nominal  3  dB  performance  loss  in  terms  of  average  -rr-  under  Rayleigh  fading  condi- 
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Figure 
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Figure  23. 
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MFSK/QPSK  Opt.  Rcvr.  Bit  Error  Performance  in  Rayleigh  Fading 
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tions,  and  performance  is  correspondingly  degraded  by  about  10  dB  at  a  BER  of  10-^  . 
For  the  M-ar\-  PSK  DBD  receivers,  performance  is  degraded  by  about  4  dB  and  by  12 
dB  at  a  BER  of  10"'  and  10-^,  respectively. 


Table  5.     EFFECT  OF  RAYLEIGH  FADING  USING  VARIOUS  RECEIVERS 

Modulation/ 
Receiver 

Rayleigh 
Fading 

-jz-  (dB)  Required  to  Achieve 
10-'  BER           10-2  BER 

2FSK  QPSK 

Optimum 
Receiver 

No 

-0.6 

03.7 

Yes 

02.9 

14.0 

2FSK  QPSK 

DBD 

Receiver 

No 

00.4 

05.6 

Yes 

03.6 

15.1 

8-PSK 

DBD 

Receiver 

No 

00.1 

03.9 

Yes 

04.4 

16.0 

16-PSK 

DBD 

Receiver 

No 

03.6 

07.4 

Yes 

07.2 

19.4 
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VI.     CONCLUSIONS 

The  results  presented  show  that  MFSK/QPSK  modulation  is  superior  to  M-ary  FSK 
in  aspects  of  error  performance  and  bandwidth  efficiency.  While  MFSK/QPSK  has  only 
slightly  better  error  performance  over  M-ary  FSK,  MFSK/QPSK  requires  only  one-third 
to  one-fourth  the  bandwidth  used  by  M-ar>'  FSK  signaling.  Also,  since  the 
MFSK/QPSK  receiver  requires  the  use  of  only  one  fourth  the  number  of  synchronized 
carriers  and  integrate  and  dump  circuits  required  by  the  M-ary  FSK  receiver,  it  can  be 
argued  that  the  MFSK/QPSK  receiver  is  less  complex  than  the  M-ary  FSK  receiver. 
The  results  also  show  that  MFSK/QPSK  is  comparable  to  M-ary  PSK  in  the  sense  that 
it  offers  improved  error  performance  over  M-ar>'  PSK,  but  does  not  reduce  bandwidth 
efficiency  as  much  as  orthogonal  signaling. 

In  comparing  the  different  types  of  MFSK,QPSK  receivers,  the  DBD  receiver  offers 
the  ability  to  decode  bits  individually  without  symbol-to-bit  mapping  circuitr}";  however, 
this  is  accomplished  at  the  expense  of  error  performance.  This  error  performance  pen- 
alty becomes  prohibitive  for  larger  values  of  k.  Since  the  advantage  of  receiver  simplicity 
is  no  longer  as  attractive  for  larger  k.  the  DBD  receiver  should  be  considered  as  an  al- 
ternative to  the  optimum  receiver  only  for  smaller  values  of  k. 

Analysis  of  the  effect  of  Rayleigh  fading  in  M-ar\'  FSK  DBD  receivers  and  in 
MFSK  QPSK  receivers  has  generally  demonstrated  severe  performance  degradation  in 
the  sense  of  producing  large  increases  in  the  bit  energy  to  noise  PSD  ratios  required  to 
maintain  only  modest  BER  performance.  The  penalties  increase  rapidly  for  lower  bit 
error  probabilities.  Unless  fading  effects  can  be  overcome,  all  of  the  receivers  studied 
become  impractical  in  the  presence  of  fading.  The  MFSK  QPSK  receivers  show  a  slight 
reduction  in  performance  degradation  over  the  M-ary  PSK  DBD  receivers.  For  the 
M-ary  PSK  receivers,  the  same  value  of  the  angle  adjustment  a  which  produced  the 
minimum  BER  for  a  non-fading  channel  also  produced  the  best  overall  performance 
under  Rayleigh  fading  conditions. 
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APPENDIX       A.   INTEGRAL  IDENTITY 

This  appendix  presents  the  proof  for  the  integral  identity  used  in  the  text  as  Eq.  4.13. 


r 


Q{<xA)Q{(iAYj,{A)dA  = 


/?tan 


-1 


J_ 

4 


'P'+     ^ 


a; 


a 


a  tan 


-1 


'     2^      1 

a   + — : 


^ 


2n   IP'  +  -^, 


In    /a^  +  -^ 


(5.1) 


We  integrate  by  parts,  i.e., 


Too  foo 


(5.2) 


with  the  following  substitutions 


u^Q{aA)Q{^A) 


{B.3) 


and 


^v  =MA)dA 
A 


Ac 


exp    - 


2Ar 


\dA 


(5.4) 


Differentiating  Eq.  B.3  with  respect  to  A,  yields 


du 
dA 


=  Q{aA)QiPA)  +  Q{oiA)Q{M) 


(5.5) 


where 
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()M)  =  -^eM) 


dA 


dA 


f*oo 


"ir^n-^r 


{B.6) 


'olA 


a 


'2n 


exp    - 


2  ^2 
O.   A 


where  the  last  equality  is  obtained  using  Liebnitz  rule  to  differentiate  the  integral.   Sim- 
ilar results  hold  for  Q[PA),  thus  yielding 


du  =  — 


0{o.A) 


A 


■'271 


exp(  -  ^j-  \  +  0{^A)  -j=-  exp( 


2  ^2 

a  A 


dA     [BJ] 


Now,  integrating  Eq.  B.4  to  obtain  v  results  in 


A 


2 


^o^n""^"" 


-exp    - 


2Ar 


(5.8) 


Solving  for  the  first  term  on  the  right  hand  side  of  Eq.  B.2  yields 


wvio   = 


-  Q{c^.A)QiPA)  exp 


A- 


2a; 


q\o) 


JL 

4 


For  the  second  term  on  the  right  hand  side  of  Eq.  B.2  ,  we  have 


(5.9) 


J     vdu=\     expf  --^ 
•'o  •'o         \      2Aq 


Q{aA) 


P 
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1271 


exp    - 


^VOT^)-=expf- 


2   .2 

a  A 


V 


'27Z 


dA 


(5.10) 


(*OOl- 


Q{y.A) 


P 


exp(  -  zl4')  +  QipA)  -^=-  exp(  -  zU') 

^!2n 


dA 
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where  2,  and  z,  are  defined  as 


U  • 


2Ar 


{B.U) 


and 


^2  = 


^+^ 


2Ar 


(5.12) 


The  form  of  these  integrals  is  similar  to  the  identity  given  by  [Ref.  9:  pp.  1-20]  as 


Erfc(a;c)  exp{b  x  )dx  = =^  In    — 


(5.13) 


where  a  and  b  may  be  complex,  and 


tan 


Real(a) 

Imag(a) 


< 


n 


(5.14) 


Using  this  identity  and  the  equality 


Q{x)  =  4-  Erfcf 


X 


V  v2 


(5.15) 


produces  an  equivalent  form,  namely 


J' 


'  2    2  1 

Q{ax)  exp{b  x  )dx  = ;=^  In 


4b^n 


a  +  ^^''2 
a  -  b^2 


(5.16) 


Using  Eq.  B.16  to  solve  Eq.  B.IO  yields 


vdu  = — 

Jq  j4nz^2 


In 


a  -y2,v2 


+ 


J4t:zJ2 


■^In 


/?  +y^2v'2_ 


(5.17) 


Focusing  on  the  argument  of  the  first  log  function  and  converting  to  complex  exponen- 
tial notation  produces 
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>^;  (X  ^li\  exp 


J  tan 


^i72 


^/a^  +  22  J  exp 


7  tan 


-^iv'2 
a 


=  exp 


=  exp 


j\  tan 


jl  tan 


-!/■  ^,72 


a 


V  ^n^2 


—  tan 


-\(  -^1^2 


a 


(5.18) 


where  the  last  equality  holds  by  restricting  both  a  and  z^  to  be  positive,  which  is  indeed 
the  case.    Using  this  result  in  Eq.  B.17  yields 


/?tan" 


\-du  = 


'^^  +  -^ 


'0 


a  tan 


-1 


In      p 


2    ,        1 


+ 


'     2^      1 


Ai 


^ 


Ar 


2n   la^+     ^ 


(5.19) 


Ai 


Finally,  substituting  Eqs.  B.9  and  B.19  into  Eq.  B.2  yields  the  right  hand  side  of  Eq.  B.l, 
thus  pro\'ing  the  identity. 
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APPENDIX       B.   MFSK/QPSK  BER  COMPUTER  PROGRAM 


* 
* 
* 
* 
* 


PROGRAM  TO  COMPUTE  THE  BIT  ERROR  RATES  FOR  THE  BI -ORTHOGONAL 
SIGNAL  SET  (BOSS)  FOR  THE  OPTIMUM  RCVR  AND  FOR  THE  MFSK/QPSK 
DIRECT  BIT  DETECTION  (DBD)  RCVR. 


PROGRAM  PRBTER 


IMPLICIT  CHARACTER'n  (A-Z) 


DOUBLE  PRECISION 
DOUBLE  PRECISION 
DOUBLE  PRECISION 
DOUBLE  PRECISION 
DOUBLE  PRECISION 
DOUBLE  PRECISION 
DOUBLE  PRECISION 
INTEGER 
INTEGER 


EBNODB,INIT, FINAL, DLTADB 

EBDVNO,SREDSI 

PRBTCR , PRBERR , INTGLl , INTGL2 , INTGL3 

DBDPBE,UBAPBE 

DC ADRE , F 1 , A 1 , B 1 , AERR , RERR , ERROR 

F2,A2,B2,F3 
PIE,SQR2PI,SQR8PI,ONE 
I,K,N,M,F,IER,NUMPTS 
LEVEL, LEVOLD 


EXTERNAL 
COMMON 


F1,F2,F3 
SREDSI,N 


INITIALIZE  CONSTANTS 


PIE 

= 

3. 141592654D0 

SQR2PI 

= 

DSQRT(2''^PIE) 

SQR8PI 

= 

DSQRT(8^'^PIE) 

K 

n 

6 

M 

= 

2-vVfK 

N 

= 

2*'HK-1) 

F 

= 

2Vr,V(K-2) 

A2 

= 

OO.ODO 

Bl 

= 

20.0D0 

B2 

= 

20.0D0 

ONE 

= 

Ol.ODO 

NUMPTS 

= 

101 

INIT 

= 

-lO.ODO 

FINAL 

= 

8.0D0 

DLTADB 

= 

(FINAL  -  INIT 

AERR 

= 

0.0 

RERR 

= 

.0000001 

LEVEL  = 

D 

•1) 


SUPPRESS  ALL  ERROR  MESSAGES. 
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c 
c 
c 


CALL  UERSETC LEVEL, LEVOLD) 
CALL  XUFLOW(O) 


COMPUTE  PROBABILITY  OF  BIT  ERROR 
FOR  VALUES  OF  EB/NO  IN  THE  RANGE  OF: 
INIT  <  EBNODB  <  FINAL 


OPTIMUM  RCVR 


BIT  ERROR 


DBD  RCVR 


BIT  ERROR 


PROBABILITY 


WRITE(8,100)  K,F 
WRITECS,'"^)'  ' 
WRITE(8,*)' 
&   UNION  BOUND' 

WRITE(8,''0' 
&   APPROXIMATION' 

VRITE(8,''0' INDEX   EB/NO  (DB)     PROBABILITY 
&   BIT  ERROR' 
EBNODB  =  INIT 
DO  10  I=1,NUMPTS 

EBDVNO  =  10.  0D0'^*(EBN0DB/10) 

SREDSI  =  DSQRT(2''^EBDVN0-"-K) 

Al     =  -SREDSI 

INTGLl  =  DCADRE(F1, A1,B1,AERR,RERR, ERROR, lER)  /  SQR2PI 

INTGL2  =  DCADRE(F2,A2,B2,AERR,RERR, ERROR, lER)  /  SQR8PI 

INTGL3  =  DCADRE(F3,A1,B2,AERR,RERR,ERR0R,IER)  /  SQR8PI 

IF( INTGL3.  GT.  0.  ODO)  THEN 

PRBTCR  =  (N-1)''-(INTGL2+INTGL3)  +  INTGLl 
ELSE 

PRBTCR  =  (N-1)*INTGL2  +  (N-1)''^INTGL3  +  INTGLl 
ENDIF 

PRBERR  =  ONE  -  PRBTCR 

DBDPBE  =  DERFC(  DSQRT(  2''^K*EBDVN0/M)  )  /  2 
UBAPBE  =  (  (M  -  2)''^DERFC(  DSQRT(K*EBDVN0/2)  ) 
&  +  DERFC(  DSQRT(K''-EBDVN0)  )   )  /4 

WRITE(8,101)  I, EBNODB, PRBERR, DBDPBE, UBAPBE 
EBNODB  =  EBNODB  +  DLTADB 
10  CONTINUE 


100  F0RMAT(2X, 'NUMBER  OF  BITS:  ', 2X, 1 1 ,6X, 'NUMBER  OF  FREQUENCIES:  ', 
&       2X,I1) 

101  F0RMAT(2X,I3,3X,F10.  5,3X,E14.  5E2,3X,E14.  5E2,3X,E14.  5E2) 
STOP 

END 


DOUBLE  PRECISION  FUNCTION  F1(Y) 


DOUBLE  PRECISION  Y, SREDSI ,SQRT2 
INTEGER  N 


COMMON 


SREDSI, N 


INITIALIZE  CONSTANT 


SQRT2  =  DSQRT(2.0D0) 
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*  COMPUTE  F1(Y,SREDSI,N) 

Fl  =  DERF(  (Y+SREDSI)/SQRT2  )''"V(n-1)  *   DEXP( -(Y**2)/2) 

RETURN 
END 


DOUBLE  PRECISION  FUNCTION  F2(Y) 
DOUBLE  PRECISION  Y,FACTR1 ,FACTR2 ,SREDSI 
DOUBLE  PRECISION  SQRT2 
INTEGER  N 

COMMON  SREDSI,N 

C  INITIALIZE  CONSTANTS 

SQRT2  =  DSQRT(2. ODO) 

C  COMPUTE  F2(Y,SREDSI,N) 

FACTRl  =  DERF(  Y/SQRT2  )'V*(N-2) 

FACTR2  =  DERF(  ( Y+SREDSI)/SQRT2  ) 

F2     =  FACTRl'VFACTR2''-DEXP(-(Y*''2)/2) 

RETURN 
END 


DOUBLE  PRECISION  FUNCTION  F3(Y) 

DOUBLE  PRECISION  Y,FACTR1 ,FACTR2 , SREDSI 

DOUBLE  PRECISION  SQRT2 

INTEGER  N 

COMMON  SREDSI, N 

C  INITIALIZE  CONSTANTS 

SQRT2  =  DSQRT(2.0D0) 

C  COMPUTE  F3(Y, SREDSI, N) 

FACTRl  =  DERF(  (Y+SREDSI)/SQRT2  )**(N-2) 

FACTR2  =  DERF(  Y/SQRT2  ) 

F3     =  FACTRl''^FACTR2'>DEXP(-(Y+SREDSI)**2/2  ) 

RETURN 
END 
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APPENDIX 


C.    16-PSK  RAYLEIGH  FADING  BER  COMPUTER 
PROGRAM 


* 

*  PROGRAM  TO  COMPUTE  THE  PROBABILITY  OF  A  BIT  ERROR  FOR  16-PSK  USING 

*  THE  DIRECT  BIT  DETECTION  (DBD)  RECEIVER  IN  THE  RAYLEIGH  FADING 

*  ENVIROMENT. 
* 

* 


PROGRAM 


DBDRAY 


IMPLICIT  CHARACTER--n  (A-Z) 


DOUBLE  PRECISION 
DOUBLE  PRECISION 
DOUBLE  PRECISION 
DOUBLE  PRECISION 
DOUBLE  PRECISION 
DOUBLE  PRECISION 
DOUBLE  PRECISION 
DOUBLE  PRECISION 
INTEGER 

EXTERNAL 
COMMON 


EBNODB , INIT , FINAL , DLTADB , PBERAY 

EBDVNO , AVESNR , ALPHA , AL$DEG , BETA 

P1BERY,V1$FTN(4) ,ARGMNT 

P2BERY,V2$FTN(4),ARG$1,ARG$2 

P3BERY, INTGRLC  8) ,TERM$1 ,TERM$2 ,TERM$3 

Fl , A , B , ERRABS , ERRREL , INTGLl , ERREST 

F2,F3,INTGL2,INTGL3 

PI,0NE,SUM$1,SUM$2 

I,INDEX,NUMPTS,IRULE 

F1,F2,F3 

AVESNR , ALPHA , BETA , INDEX 


CALL  ERSET(3,0,-1) 


PI 

ONE 

AL$DEG  = 
ALPHA  = 
BETA 

V1$FTN(1) 
V1$FTN(2) 
V1$FTN(3) 
V1$FTN(4) 


3. 141592654D0 

l.ODO 

14. 25D0 

AL$DEG  *   PI/ 180 

PI/4  -  ALPHA 

=  DCOS( ALPHA) 
=  DCOS(BETA) 
=  DSIN( ALPHA) 
=  DSIN(BETA) 


SUPPRESS  IMSLIO  ERROR  MESSAGES. 


INITIALIZE  CONSTANTS 


V2$FTN(1)  =  DCOS( ALPHA)  -  DSIN( ALPHA) 
V2$FTN(2)  =  DCOS(ALPHA)  +  DSIN(ALPHA) 
V2$FTN(3)  =  DCOS(BETA)  -  DSIN(BETA) 
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V2$FTN(4)  =  DCOS(BETA)  +  DSIN(BETA) 


NUMPTS 
INIT 
FINAL 
DLTADB 

A 

B 

ERRABS 
ERRREL 
IRULE 


101 

■10.  ODO 
20.  ODO 
(FINAL 

O.ODO 
PI/8 
O.ODO 
0.  0001 
2 


-  INIT) /( NUMPTS- 1) 


C  COMPUTE  PROBABILITY  OF  BIT  ERROR 

C  FOR  VALUES  OF  EBNODB  IN  THE  RANGE  OF: 

C  INIT  <  EBNODB  <  FINAL 

VRITE( 8,*) 'EFFECT  OF  RAYLEIGH  FADING  ON  16-PSK,  DBD  RCVR' 

WRITE(8,100)  AL$DEG 

WRITE(8,*)'  ' 

WRITE (8,'^)'  DBD  RCVR' 

WRITE(8,''-)'  BIT  ERROR' 

VRITE(8,'^)'INDEX  EB/NO  (DB)     PROBABILITY' 

EBNODB  =  INIT 


BEGIN  MAIN  LOOP. 
DO  70  1=1, NUMPTS 
EBDVNO  =  10.0D0--'^'^(EBN0DB/10) 
AVESNR  =  2  'V  EBDVNO 

CALCULATE  PROBABILITY  OF  IB  ERROR 
SUM$1  =  O.ODO 

DO  10  INDEX=1,4 
ARGMNT  =  EBDVNO  /  (4  *   EBDVNO  *   V1$FTN(  INDEX) *-'^2  +  ONE) 
SUM$1   =  SUM$1   +  V1$FTN( INDEX)  *   DSQRT( ARGMNT) 
10    CONTINUE 

PIBERY  =  (2. ODO  -  SUM$l)/4 

CALCULATE  PROBABILITY  OF  2B  ERROR 
SUM$1  =  O.ODO 

DO  20  INDEX=1,4 
ARG$1  =  EBDVNO  /  (EBDVNO  *  V2$FTN( INDEX)**2  +  ONE/2) 
SUM$1  =  SUM$1  +  V2$FTN(  INDEX)  '>  DSQRT(ARG$1) 
20    CONTINUE 

SUM$2  =  0.  ODO 

DO  30  INDEX=1,4 
ARG$1  =  DSQRT(   V2$FTN( INDEX)**2  +  ONE  /  (2*EBDVN0)   ) 
ARG$2  =  ARG$1  /  V2$FTN(  INDEX  -  ( -1)*^( INDEX)  ) 
SUM$2  =  SUM$2  +  V2$FTN( INDEX)  *  DATAN(ARG$2)  /  ARG$1 
30    CONTINUE 

P2BERY  =  ONE/2  -  SUM$l/4  +  SUM$2/(2*PI) 

CALCULATE  PROBABILITY  OF  3B  ERROR 
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DO  40   INDEX=1,8 
CALL  DQDAG   (Fl ,A,B,ERRABS,ERRREL, IRULE, INTGL1,ERREST) 
TERM$1   =  DSQRT(AVESNR)'VINTGLl/4 

CALL  DQDAG   (F2,A,B ,ERRABS ,ERRREL,IRULE, INTGL2,ERREST) 
TERM$2  =   INTGL2 

CALL  DQDAG   (F3,A,B,ERRABS ,ERRREL, IRULE, INTGL3,ERREST) 
TERM$3   =  DSQRT(AVESNR)''-INTGL3 

INTGRL( INDEX)   =  TERM$1    -    (TERM$2  +  TERM§3)/(4*PI) 
40  CONTINUE 

SUM$1  =  O.ODO 

DO  50   INDEX=5,8 
SUM$1  =  SUM$1  +   INTGRLC INDEX) 
50  CONTINUE 

SUM$2  =  0.  ODO 

DO   60   INDEX=1,4 
SUM$2  =  SUM$2  +   INTGRLC INDEX) 
60  CONTINUE 

P3BERY  =  ONE/2  +  SUM$1    -   SUM$2 
PBERAY  =  PlBERY/2  +  P2BERY/4  +  P3BERY/4 
WRITE(8,101)    I, EBNODB, PBERAY 
EBNODB   =  EBNODB   +  DLTADB 
70  CONTINUE 

C  END  MAIN  LOOP 

100  F0RMAT(2X, 'ALPHA  =' , 1X,F5. 2, IX, 'DEGREES. ' ) 

101  F0RMAT(2X,I3,3X,F10.  5,3X,E14.  5E2,3X,E14.  5E2,3X,E14.  5E2) 

STOP 

END 


4\    *\    *\    4%    4\    4\    4\    4\    «*  4%    4\    4\    4\    4\    4\    *\    4\    4\    4\    4\    4\    4\    t\    4\    4\    4\    4\    4\    4%    4\    4\    4\    4\    4\    4\    4\    4\    #h  **  4\    4\    4\    4\    4\4\    4\    4\7\4T  4\    4\    4\    4%    4%    4^    4\    4\    *S  *V  4i  4\    4\    4\    4\f\    4\    4\    *%    4\4\    4\    *f 

DOUBLE   PRECISION  FUNCTION     Fl( GAMMA) 

IMPLICIT  CHARACTER^n    (A-Z) 
DOUBLE   PRECISION     GAMMA, AVE SNR, ALPHA, BETA 
DOUBLE   PRECISION     TERM$1 ,FINDEX,GINDEX 
INTEGER  INDEX 

COMMON  AVE SNR, ALPHA, BETA, INDEX 

C  COMPUTE  Fl( INDEX, GAMMA, AVESNR, ALPHA, BETA) 

C 

GO  TO   (10,20,30,40,50,60,70,80)    INDEX 

10  FINDEX  =  DSIN( GAMMA  -  ALPHA) 
GINDEX  =  DCOS( GAMMA  -  ALPHA) 
GO  TO   100 
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20  FINDEX  =  DCOS( GAMMA  -  ALPHA) 
GINDEX  =  DSIN( GAMMA  -  ALPHA) 
GO  TO   100 

30  FINDEX  =  DCOS( GAMMA  +  ALPHA) 
GINDEX  =  DSIN( GAMMA  +  ALPHA) 
GO  TO   100 

40  FINDEX  =  DSIN( GAMMA  +  ALPHA) 
GINDEX  =  DCOS( GAMMA  +  ALPHA) 
GO  TO   100 

50  FINDEX  =  DSIN( GAMMA  -  BETA) 
GINDEX  =  DCOS( GAMMA  -  BETA) 
GO  TO   100 

60  FINDEX  =  DCOS( GAMMA  -  BETA) 
GINDEX  =  DSIN( GAMMA  -  BETA) 
GO  TO   100 

70  FINDEX  =  DCOS( GAMMA  +  BETA) 
GINDEX  =  DSIN( GAMMA  +  BETA) 
GO  TO   100 

80  FINDEX  =  DSIN( GAMMA  +  BETA)     - 
GINDEX  =  DCOS( GAMMA  +  BETA) 

100  TERM$1  =  AVESNR  ''^  FINDEX^'-"2        +      1.  ODO/2 
Fl  =  GINDEX  /   TERM$1   /   DSQRT(TERM$1) 

RETURN 

END 


<J^*J*  y*  ^^  ^^  «'-  -f'  ^-  -';•  >'^  -*^  o'-  -^^  <^^  -''  •^-  »'-  -^-  o^--  •^-  -f'  •;''  ^-  •-*'  -''  •^'■-  -*'  -*-•  •J'  >'-  ■^*«  <J-  -J'  J*  •!'  J'  -i^  -*-  -''  »%  J*  -f'  -J-  *''  -J'  -'c  -I-  •J'  -*'  -Jf  -J'  -''■  oL.  J.  .J^  ^^  kt«  ^«.t^  JU  J^  -J'  -f-  <J'  ^«  ^^  JU  JU  y^  .J^  y^  ^ 

DOUBLE   PRECISION  FUNCTION     F2(GAMMA) 

IMPLICIT  CHARACTER-a   (A-Z) 
DOUBLE   PRECISION     GAMMA, AVESNR, ALPHA, BETA 
DOUBLE   PRECISION     TERM$1,TERM$2 , FINDEX, GINDEX 
INTEGER  INDEX 

COMMON  AVESNR, ALPHA, BETA, INDEX 

C  COMPUTE  F2 (INDEX, GAMMA, AVESNR, ALPHA, BETA) 

C 

GO  TO   (10,20,30,40,50,60,70,80)    INDEX 

10  FINDEX  =  DSIN( GAMMA  -  ALPHA) 
GINDEX  =  DCOS( GAMMA  -  ALPHA) 
GO  TO   100 

20  FINDEX  =  DCOS( GAMMA   -   ALPHA) 
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GINDEX  =  DSINC GAMMA   -  ALPHA) 
GO  TO   100 

30  FINDEX  =  DCOS( GAMMA  +  ALPHA) 
GINDEX  =  DSIN( GAMMA  +  ALPHA) 
GO  TO   100 

40  FINDEX  =  DSIN( GAMMA  +  ALPHA) 
GINDEX  =  DCOS( GAMMA  +  ALPHA) 
GO  TO   100 

50  FINDEX  =  DSINC GAMMA  -  BETA) 
GINDEX  =  DCOS( GAMMA  -  BETA) 
GO  TO   100 

60  FINDEX  =  DCOS( GAMMA  -  BETA) 
GINDEX  =  DSINC GAMMA  -  BETA) 
GO  TO   100 

70  FINDEX  =  DCOSC GAMMA  +  BETA) 
GINDEX  =  DSINC GAMMA  +  BETA) 
GO  TO   100 

80  FINDEX  =  DSINC GAMMA  +  BETA) 
GINDEX  =  DCOSC GAMMA  +  BETA) 

100   TERM$1   =  AVESNR  ^>-   FINDEX^'^"-2        +      1.  ODO/2 
TERM$2   =  AVESNR  '''   GINDEX-----2 
F2  =  TERM$2   /    CTERM$1---TERM$2  +  TERM$1*"2) 

RETURN 
END 


t\     *%     *\    0\    4\     ««     «*     #t    *\     «»     4%     t\     *\     *\     f\     4\     t\     4\    ¥\     »\    t\     t\     t\     t\     4\     *y     4\     *\     *\     4\     4\     »\     t\     *\     4%    4\     4\     4\     4\4Z4\     4\7\    ^7  4^7%   4%    4S     4\  7Z  4*  TT  4\  TZ  4\     4\4V4%4%7V4*     4\     4k     4h     4\     4\     4%    4%     4\    4\     4\  TV 

DOUBLE   PRECISION  FUNCTION     F3C GAMMA) 

IMPLICIT  CHARACTER-'-l    CA-Z) 

DOUBLE   PRECISION      GAMMA, AVESNR, ALPHA, BETA 

DOUBLE   PRECISION     TERM$1 ,TERM$2 ,FACTR1 ,FACTR2 , FINDEX, GINDEX 

INTEGER  INDEX 

COMMON  AVESNR, ALPHA, BETA, INDEX 

C  COMPUTE  F3C INDEX, GAMMA, AVESNR, ALPHA, BETA) 

C 

GO  TO   CIO, 20, 30, 40, 50, 60, 70, 80)    INDEX 

10  FINDEX  =  DSINC GAMMA  -  ALPHA) 
GINDEX  =  DCOSC GAMMA  "  ALPHA) 
GO  TO    100 

20  FINDEX  =  DCOSC GAMMA   -   ALPHA) 
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GINDEX  =  DSINC GAMMA  -  ALPHA) 
GO  TO  100 

30  FINDEX  =  DCOS( GAMMA  +  ALPHA) 
GINDEX  =  DS IN (GAMMA  +  ALPHA) 
GO  TO  100 

40  FINDEX  =  DSIN( GAMMA  +  ALPHA) 
GINDEX  =  DCOS( GAMMA  +  ALPHA) 
GO  TO  100 

50  FINDEX  =  DSIN( GAMMA  -  BETA) 
GINDEX  =  DCOS( GAMMA  -  BETA) 
GO  TO  100 

60  FINDEX  =  DCOS( GAMMA  -  BETA) 
GINDEX  =  DSIN( GAMMA  -  BETA) 
GO  TO  100 

70  FINDEX  =  DCOS( GAMMA  +  BETA) 
GINDEX  =  DSIN( GAMMA  +  BETA) 
GO  TO  100 

80  FINDEX  =  DSIN( GAMMA  +  BETA) 
GINDEX  =  DCOS( GAMMA  +  BETA) 

100  TERMSl  =  AVESNR  *   FINDEX''^'V2   +   1.  ODO/2 
TERM$2  =  AVESNR  '"^  GINDEX''^''-2 
FACTRl  =  GINDEX  /  TERM$1  /  DSQRT(TERM$1) 
FACTR2  =  DATAN(  DSQRT(TERM$2/TERM$1)  ) 
F3     =  FACTRl  ■>'■   FACTR2 


RETURN 
END 
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APPENDIX       D.   MFSK/QPSK  RAYLEIGH  FADING  BER  COMPUTER 

PROGRAM 


* 

*  PROGRAM  TO  COMPUTE  THE  BIT  ERROR  RATES  FOR  THE  B I -ORTHOGONAL 

*  SIGNAL  SET  (BOSS)  FOR  THE  OPTIMUM  RCVR  WHEN  SUBJECTED  TO  RAYLEIGH 

*  FADING. 
* 


PROGRAM  PBTERA 


IMPLICIT  CHARACTER^-1  (A-Z) 


DOUBLE  PRECISION 
DOUBLE  PRECISION 
DOUBLE  PRECISION 
DOUBLE  PRECISION 
DOUBLE  PRECISION 
INTEGER 
INTEGER 

EXTERNAL 
COMMON 


EBNODB , INIT, FINAL, DLTADB 

EBDVNO,SOSQRD,S 

PBCRAY,PBERAY 

DCADR2 , INTGND , LIL , UIL , AERR , RERR , ERROR 

ONE 

I,K,N,M,F,IER,NUMPTS 

LEVEL, LEVOLD 

INTGND 
SOSQRD,S,N 


INITIALIZE  CONSTANTS 


K 

= 

5 

M 

= 

2V«VK 

N 

= 

2Vryc(K-l) 

F 

= 

2V,-Vr(K-2) 

LIL 

= 

00.  ODO 

ONE 

= 

Ol.ODO 

NUMPTS 

= 

101 

INIT 

= 

-10.  ODO 

FINAL 

= 

20.  ODO 

DLTADB 

= 

(FINAL  - 

AERR 

= 

0.  0 

RERR 

= 

.0001 

INIT) /( NUMPTS- 1) 


SUPPRESS  ALL  ERROR  MESSAGES. 


LEVEL  =  0 

CALL  UERSET( LEVEL, LEVOLD) 

CALL  XUFLOW(O) 
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c 
c 
c 


COMPUTE  PROBABILITY  OF  BIT  ERROR 
FOR  VALUES  OF  EB/NO  IN  THE  RANGE  OF: 
INIT  <  EBNODB  <  FINAL 


WRITE(8 
WRITE (8 
WRITE(8 
WRITE (8 
WRITE(8 
WRITE(8 
WRITE (8 
EBNODB  = 


*) 'EFFECT  OF  RAYLEIGH  FADING  ON  BOSS' 

100)  K,F 
*)  ' 

"0' INDEX 
INIT 


EB/NO  (DB) 


OPTIMUM  RCVR' 
BIT  ERROR 
PROBABILITY  ' 


DO  10  I=1,NUMPTS 

EBDVNO  =  10.0D0**(EBN0DB/10) 
SOSQRD  =  EBDVNO  *  K 

UIL    =  MIN(18.5D0,10*DSQRT( SOSQRD)  ) 

PBERAY  =  DCADR2( INTGND, LIL, UIL, AERR,RERR, ERROR, lER)  /  SOSQRD 
WRITE(8,101)  I, EBNODB, PBERAY 
EBNODB  =  EBNODB  +  DLTADB 
10  CONTINUE 

100  F0RMAT(2X, 'NUMBER  OF  BITS:  ' ,2X, II ,6X, 'NUMBER  OF  FREQUENCIES:  ', 
&       2X,I1) 

101  FORMAT(2X,I3,3X,F10.5,3X,E14. 5E2,3X,E14.  5E2,3X,E14.  5E2) 

STOP 
END 


DOUBLE  PRECISION  FUNCTION  INTGND( SREDSI) 
IMPLICIT  CHARACTER^'- 1  (A-Z) 


DOUBLE  PRECISION 
DOUBLE  PRECISION 
DOUBLE  PRECISION 
DOUBLE  PRECISION 
DOUBLE  PRECISION 


EBDVNO , SREDSI , SOSQRD , S , THRHLD 
PRBTCR , PRBTER , ORTERM , APTERM 
INTGL1,INTGL2,INTGL3 
DCADRE,F1,A1,B1,AERR,RERR,ERR0R 
F2,A2,B2,F3 


DOUBLE  PRECISION  PIE,SQR2PI ,SQR8PI ,ONE 


INTEGER 

EXTERNAL 
COMMON 


N,IER 

F1,F2,F3 
SOSQRD, S,N 


INITIALIZE  CONSTANTS 


THRHLD 

PIE 

ONE 

SQR2PI 

SQR8PI 

AERR 

RERR 


6.2D0 

3. 141592654D0 

l.OODO 

DSQRT(2*PIE) 

DSQRT(8^-PIE) 

0.  0 

. 0000001 
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A2 

=  00.  ODO 

Bl 

=   20.  ODO 

B2 

=   20.  ODO 

c 
c 

COMPUTE  PROBABILITY  OF  BIT  ERROR  AS  A 
FUNCTION  OF  SREDSI  AND  N 
IF(SREDSI.GT.  THRHLD)  THEN 

USE  UNION  BOUND  APPROXIMATION 
ORTERM  =  (2*N  -  2)'''DERFC(  SREDSI/2  ) 
APTERM  =  DERFCC  SREDSI  /  DSQRT(2.0D0)  ) 
PRBTER  =  (  ORTERM  +  APTERM  )  /  4 
ELSE 

USE  EXACT  CALCULATION 
S      =  SREDSI 
Al     =  -SREDSI 

INTGLl  =  DCADRE(F1,A1,B1,AERR,RERR,ERR0R,IER)  /  SQR2PI 
INTGL2  =  DCADRE(F2,A2,B2,AERR,RERR,ERR0R,IER)  /  SQR8PI 
INTGL3  =  DCADRE(F3,A1,B2,AERR,RERR, ERROR, lER)  /  SQR8PI 
IF( INTGL3. GT. 0. ODO)  THEN 

PRBTCR  =  INTGLl  +  (N-l)*( INTGL2+INTGL3) 
ELSE 

PRBTCR  =  INTGLl  +  (N-1)'''INTGL2  +  (N-1)*INTGL3 
ENDIF 

PRBTER  =  ONE  -  PRBTCR 

ENDIF 

COMPUTE  INTGNDC SREDSI, S0SQRD,N). 

INTGND  =  PRBTER  *  SREDSI  '"^  DEXP(  -(SREDSI*'>2)  /  (2  *  SOSQRD)  ) 

RETURN 
END 


DOUBLE  PRECISION  FUNCTION  F1(Y) 

DOUBLE  PRECISION  Y, SOSQRD, S,SQRT2 

INTEGER  N 

COMMON  SOSQRD, S,N 

C  INITIALIZE  CONSTANT 

SQRT2  =  DSQRT(2. ODO) 

C  COMPUTE  F1(Y,S,N) 

Fl  =  DERF(  (Y+S)/SQRT2  )**(N-1)  *  DEXP( -(Y**2)/2) 

RETURN 
END 
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DOUBLE  PRECISION  FUNCTION  F2(Y) 

DOUBLE  PRECISION  Y,FACTR1 ,FACTR2 ,SOSQRD, S 

DOUBLE  PRECISION  SQRT2 

INTEGER  N 

COMMON  SOSQRD,S,N 

C  INITIALIZE  CONSTANTS 

SQRT2  =  DSQRT(2.  ODO) 

C  COMPUTE  F2(Y,S,N) 

FACTRl  =  DERF(  Y/SQRT2  )**(N-2) 

FACTR2  =  DERF(  (Y+S)/SQRT2  ) 

F2     =  FACTRl'>FACTR2''-DEXP(-(Y**2)/2) 

RETURN 

END 


DOUBLE  PRECISION  FUNCTION  F3(Y) 

DOUBLE  PRECISION  Y,FACTR1 ,FACTR2,S0SQRD,S 

DOUBLE  PRECISION  SQRT2 

INTEGER  N 

COMMON  SOSQRD,S,N 

C  INITIALIZE  CONSTANTS 

SQRT2  =  DSQRT(2.0D0) 

C  COMPUTE  F3(Y,S,N) 

FACTRl  =  DERF(  (Y+S)/SQRT2  )-'"V(N-2) 

FACTR2  =  DERF(  Y/SQRT2  ) 

F3     =  FACTRl'>FACTR2*DEXP(-(Y+S)**2/2  ) 

RETURN 
END 
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